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COMPARISONS  WITH  BALLOON  DATA 
by 

Hsiao-hua  Kuo  Burke 


ABSTRACT 


Measurements  of  air  conductivity  have  indicated  that  the  conduc- 
tivity of  the  atmosphere  can  be  approximated  as  an  exponentially  increas- 
ing function  of  the  altitude  with  a scale  height  of  5-6  km.  Because  of 
this  property,  the  atmospheric  electric  fields  produced  by  a variety  of 
charge  sources  cannot  be  described  by  the  simple  Coulomb  field  relation- 
ship. In  this  thesis,  a new  theoretical  approach  is  employed  to  calcu- 
late the  large-scale  weather-induced  atmospheric  fields  at  different 
altitudes.  The  Green's  function  method  is  used  to  obtain  analytic 
solutions  for  different  charge  sources.  This  mapping  technique  is  then 
extended  into  the  lower  ionosphere  where  the  conductivity  is  of  a tensor 
form. 

The  model  is  used  to  evaluate  high-altitude  electric-field  data 
from  a balloon  flight.  The  balloon  was  launched  on  August  3,  1969,  from 
Penhold,  Canada,  and  measured  vector  electric  fields  at  28  km  for  a 
period  of  24  hours.  Data  included  the  electric  fields  of  a local  thun- 
derstorm during  the  earlier  portion  of  the  flight  and  fair  weather 
fields  during  the  later  portion. 

Comparisons  between  theoretical  models  and  balloon  data  indicate 
that  during  the  storm  the  balloon  was  about  22  km  south  of  the  thunder- 
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CHAPTER  1 


INTRODUCTION 

I . General  Review  of  Atmospheric  Electricity 

For  a long  time,  the  study  of  atmospheric  electricity  was  con- 
fined to  the  theory  of  "concentric  spherical  capacitor"  (cf.  Chalmers, 
1967).  The  lower  boundary  of  this  capacitor  was  the  earth's  surface 
and  the  electrosphere,  a highly  conducting  layer  at  about  50-70  km,  was 
the  upper  boundary.  The  electrosphere  was  defined  as  the  height  at  which 
the  equipotential  condition  was  reached.  Being  a good  conductor,  it 
acted  as  a perfect  electrostatic  shield.  As  a result,  thunderstorms  and 
other  electrical  phenomena  in  the  lower  atmosphere  would  not  have  any 
effect  outside  the  electrosphere.  At  the  same  time,  charges  bombarding 
the  ionosphere  from  outside  should  not  have  any  effect  in  the  lower 
atmosphere  unless  they  penetrate  through  the  electrosphere.  In  this 
classical  picture,  the  thunderstorm  activity  would  act  as  the  major 
mechanism  for  charge  generation  and  the  fair  weather  conduction  current 
as  the  consumer. 

This  "close-in"  theory  on  atmospheric  electricity  works  in  most 
aspects  except  it  does  not  allow  any  horizontal  atmospheric  electric 
field  which  is  known  to  exist  both  from  weather  and  ionospheric  origins. 

A closer  look  at  the  electrosphere  also  indicates  that  such  an  ideal 
equipotential  layer  does  not  exist.  Thus  the  perfect  capacitor  now 
becomes  a leaky  one  which  allows  the  atmospheric  electric  phenomena  to 
interact  with  the  ionosphere  and  magnetosphere.  Influences  of  solar  and 
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the  air  so  that 
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a ~ n+ek+  + n-ek. 


(1-2) 


f 
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i 

1 
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where  n+  and  n_  are  number  densities  of  positive  and  negative  small 

ions  and  k+  and  k_  are  their  mobilities.  Near  earth's  surface, 

n+  = 600  cm  3 , n_  ~ 500  cm-3,  k+  ~ 1.3  cm2V-1  sec-1  and  k_  ~ 1.6  cm2V-1 

sec-1  which  give  a conductivity  of  about  2.5  x 10-1 ‘'ohm-^-1  or 

2.5  x 10~3 &olun-1cra-1 . 

The  conductivity  of  air  increases  rapidly  with  altitude  due  to 
the  following  reasons: 

(1)  The  mobility  of  ions  increases  when  the  neutral  number  density 
decreases . This  is  the  dominant  factor  for  the  increase  of 
conductivity.  To  prove  this  property,  we  start  from  the  relation 
for  the  ion  mobility  k. 


k 


where  the  parameters  are 

v,  the  mean  thermal  velocity  of  the  ions, 

1,  the  mean  free  path  of  ions,  and 
T,  the  temperature  of  the  air. 

The  increase  of  conductivity  with  respect  to  v and  1 is 
obvious.  Conductivity  also  increases  with  the  decrease  of 
viscosity  and  in  a gas  viscosity  increase  proportionally  with 
temperature  which  implies  the  relationship  of  k ® ^ . Since 
1 a — , n being  the  neutral  number  density,  and  v /f,  we 


get  k 


OC 


_1_ 

n^T 


. Then  from 
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-z/h 

n “ e 

where  z is  the  altitude  and  h the  scale  height  of  density, 
we  reach  the  conclusion  that  conductivity  increases  as  we  go 
higher.  Figure  1-1  shows  the  variation  of  ion  mobility  k 
with  atmospheric  neutral  number  density. 

(2)  The  presence  of  numerous  large  aerosols  lowers  the  conductivity 
of  air  near  ground.  After  small  ions  attach  themselves  to  aero- 
sols, their  mobility,  and  hence  the  conductivity  of  air,  is 
decreased  due  to  the  increase  of  the  collisional  cross  sectional 
area. 

(3)  The  effect  of  cosmic  radiation  on  ionozation  increases  with  alti- 
tude (Table  1-1)  in  the  lower  atmosphere. 


III.  The  Conductivity  Profile 

(i)  Atmospheric  Conductivity  Profile 


Variations  of  electrical  conductivity  in  the  atmosphere  as  a func- 
tion of  altitude  have  been  observed  by  a number  of  people  (cf.  Benbrook 
et  al.  1974;  Bragin,  1967;  Paltridge,  1965;  Hatakeyama,  1963;  Kraakevik, 
1958) . Cole  and  Pierce  (1965)  deduced  a conductivity  profile  versus 
height  (Figure  1-2)  based  on  observations  on  variations  of  electron  and 
ion  concentrations  and  effective  collisional  frequencies  with  respect  to 
altitude.  From  this  figure  the  conductivity  profile  in  the  atmosphere 
was  modeled  as 
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o ( z) 


3 x 10  1 5 exp(z/0.82)  ohm  1m_1 

0 < z < 3.6 

1 x 10-13  exp(z/4.1)  ohm-1m_1 

3.6  < z £ 17.7 

6 x 10-13  exp(z/7.0)  ohm-1m-1 

17 . 7 < z < 40 


where  z is  the  height  in  kilometers. 

Because  of  this  variational  property  of  conductivity,  the  air 
can  not  be  treated  as  a uniformly  conducting  medium  when  the  electric 
field  of  interest  is  extended  beyond  a few  kilometers. 

As  a convenient  large-scale  working  model,  the  conductivity 
profile  can  be  approximated  as 


a(z)  = aQ  exp(2kz)  (1-3) 

where  (2k)-1  is  called  the  scale  height  of  conductivity.  The  accepted 
values  of  the  scale  height  vary  between  5 and  6 km. 

Table  1-2  shows  conductivities  at  different  altitudes  calculated 
from  (1-3)  with  different  scale  heights.  Figure  1-3  shows  the  computed 
values  as  compared  to  the  Cole  and  Pierce  profile  and  other  measurements 


(ii)  Ionospheric  Conductivity  Profile 

The  region  above  70  km  is  called  the  ionosphere.  The  ionization 
sources  in  this  region  are  different  from  those  in  the  atmosphere.  During 
daytime,  the  major  source  of  ionization  is  the  flux  of  the  extremely  high 


Alt. 

4 (ohm 

-1  -U 
m ) 

(km) 

| (2k)_1=  5 km 

( 2k)  6 km 

0 

>0x  10'14 

2. 2x  10 


1.6x  10" 


8.4x10" 


4.4x10" 


Table  1-2 

Computed  conductivities  of  air  at  different 
altitudes  based  on  (l-3) 


I 


energy  solar  photons.  During  nighttime,  the  major  source  is  the  precipi- 
tation of  magnetospheric  particles.  In  addition,  since  the  electron  gyro- 
frequency  becomes  comparable  to  the  collisional  frequency,  magnetic  forces 
have  to  be  taken  into  consideration.  The  conductivity  of  a plasma  that 
has  magnetic  field  lines  imbedded  in  it  is  no  longer  isotropic;  it  is  a 
tensor  rather  than  a scalar. 


o=o  a,  0 
: 2 1 


(1-4) 


(choosing  z to  point  along  the  direction  of  B) . 


As  for  a two  component  plasma  (electrons  and  one  species  of  posi- 


tive ion)  : 


eo“e  Ve 

°1  " 2 2 
ve  ^eH 


£ow+  v+ 

VHh 


is  called  the  "Pederson  conductivity;" 


a = - + £oa)+2a)+H 

v 2 -hi)2  v 2-hu2 
e eH  + eH 


is  called  the  "Hall  conductivity;"  and 


2 2 

eo“e  EoV 

ao  = + 

Ve  V+ 


is  the  "specific  conductivity"  where 
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<jd,u  = e B/m,  is  the  gyrof requency  of  ions, 

"rrl  + 

W+  = Nke2/m+eo  is  the  plasma  frequency  of  ions, 

v+  is  the  collision  frequency  for  ions 

and  similarly  for  w u)  and  v of  electrons. 

eH  e e 

Calculated  values  of  a , 0,  and  a at  different  altitudes  in 

0 i 2 

the  ionosphere  can  be  found  in  detail  in  the  Satellite  Environment  Hand- 
book (1965).  The  necessary  approximations  for  field  modeling  will  be 
discussed  in  Chapter  3.  Figure  1-4  is  a sample  conductivity  profile  up 
to  ^150  km.  Above  150  km,  the  specific  conductivity  becomes  so  large 
that  magnetic  field  lines  can  be  considered  as  equipotentials . 


Altitude  (km) 
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Figure  1-2 


Profiles  of  conductivity  versus  altitude 


Suffixes  +,  -,  and  e refer  to  positive  ions, 

negative  ions,  and  electrons  respectively.  (Cole  and  Pierce,  1965) 
In  this  Figure,  6^,  6 and  A + are  obtained  by  the  measured 
values  of  electron  and  ion  concentrations,  their  effective 
collisional  frequencies  and  the  relation 

<4  = n e /m  11 

where  n,  e,  m,i/are  number  density,  electric  charge,  mass 
and  collisional  frequency  respectively.  A total  ^en 
the  summation  of  A , A+  and  A • 


I I >. 'H  .1...  i i .iwnui  i; 
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Figure  1-3 

Computed  values  of  conductivities  as  compared  with  the 

Cole  and  Pierce  profile  (I965)  and  other  measurements  by 

Mozer  and  Serlin  (1969)  and  Paltridge  (1965)* 

The  conductivity  of  air  near  ground  by  Cole  and  Pierce 

profile  is  higher  than  the  average  accepted  value 
_14  _i  _i 

( 3x10  ohm  m ) due  to  the  neglegence  of  the  existence 
of  large  ions,  which  are  formed  as  a result  of  local 
contamination  of  the  atmosphere  and  reduce  the  conductivity 
near  the  earth's  surface. 
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Figure  1-4 

A representative  conductivity  profile  with  an  isotropic 
(scalar  conductivity)  region  below  70  km  and  an 
anisotropic  (tensor  conductivity)  region  above. 

(Park  and  Dejnakarintra,  1973) 
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CHAPTER  2 


THEORETICAL  MODEL  OF  FIELD  MAPPING 
IN  THE  ATMOSPHERE 


I . Summary  of  Work  in  the  Past 

The  purpose  of  this  modeling  of  atmospheric  electric  field  mapping 
to  different  altitudes  is  three-fold.  The  general  formulation  for  field 
mapping  is  aimed  such  that: 

(1)  it  can  apply  not  only  to  the  atmosphere  but  also  to  the  ionosphere, 

(2)  it  can  be  used  for  both  fair  weather  and  disturbed  weather  conditions, 
and 

(3)  it  should  be  possible  to  consider  different  geometries  of  charge 
sources  for  thunderstomr  field  mapping.  For  example,  in  a thunder- 
cloud model,  simple  point  charge  sources  (monopole,  dipole,  etc.) 
are  sufficient  for  examining  field  variations  at  large  distances. 
However,  disc,  cylinder  or  other  shapes  of  charge  sources  need  to 
be  used  in  order  to  study  neary-by  field  changes. 

The  research  of  large  scale  mapping  of  atmospheric  electric  field 
that  considered  the  variational  property  of  air  conductivity,  i.e.,  the 
conductivity  increases  with  altitude  in  an  exponential  manner,  can  be 
summarized  by  the  following  representative  papers: 

(1)  Holzer  and  Saxon  (1952)  and  Anderson  and  Freier  (1969).  Holzer  and 
Saxon  were  the  first  ones  who  employed  the  property  of  exponentially 
increasing  conductivity  of  air  to  calculate  the  electric  field  and 
conduction  currents  at  different  distances  due  to  a simpler  dipolar 
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thundercloud.  Anderson  and  Freier  added  in  possible  time  variation 
on  potential  changes.  However,  their  calculations  were  limited  to 
the  atmosphere;  in  other  words,  the  range  of  validity  does  not  extend 
into  the  ionosphere  where  the  conductivity  is  of  a tensor  form  rather 


than  a scalar. 

(2)  Park  and  Dejnakarintra  (1973).  The  mapping  of  thundercloud  electric 


field  at  middle  and  higher  latitudes  was  investigated  as  a three- 
dimensional  boundary  problem.  The  electrical  conductivity  was  repre- 
sented by  several  piecewise  exponential  functions  of  altitude,  and 


the  anistropy  of  air  was  taken  into  account  above  70  km.  The  geo- 
magnetic field  lines  were  assumed  to  be  straight  and  vertical  below 
150  km  in  the  regions  considered.  They  derived  solutions  for  both 


the  atmosphere  and  ionosphere  but  did  not  achieve  an  analytic  solution. 
As  a result,  involved  numerical  integrations  had  to  be  used.  Further- 


more, the  results  were  good  only  for  point  charge  sources. 


II.  Basic  Equations 


Assuming  an  exponential  increase  in  conductivity  with  respect  to 
altitude,  the  basic  equations  governing  electric  field  in  the  atmosphere 
are  the  following  ( (2— 1 ) — ( 2—5) ) : 


a = O exp(2kz) 
z 0 


(2-1) 


x e = 0 or  E = - V4> 


(2-2) 


where  0 and  0 are  conductivities  of  air  at  altitude  z and  ground 

e 0 
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respectively,  (2k)  is  the  scale  height  of  conductivity,  E is  the 

electric  field  intensity  in  the  air  and  $ is  the  electric  potential 

with  respect  to  the  earth  surface. 

In  steady  state  conditions  V • i ^ , = 0 where  i , = i +i  , 
J J total  J total  Ja  Js’ 

the  summation  of  current  density  in  the  air  and  that  at  the  source  region 

due  to  convection.  We  then  have 


j = OE 
a 


and 


• j = 0 or  V • j = -V  • i 

total  Ja  s 


(2-3) 


(2-4) 


And  Poisson's  equation  is 


V • E = P/fe0  or  V24>  = _p>fe0 


(2-5) 


where  p is  the  space  charge  density  and  is  the  permittivity  of  air. 

Cor  bination  of  (2-1),  (2-2),  (2-3),  (2-4)  gives 


Vz  + 2k 


<t>  = $ • j 


(2-6) 


Substituting  (2-5)  into  (2-6) , we  obtain 


The  first  term  is  space  charge  in  the  air 


(2-7) 


p . 

air 


2ke 


3$ 
0 9z 


(2-8) 


The  second  term  is  zero  everywhere  else  except  at  the  source  of  charge 
such  that 
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$ • V 


(2-9) 


From  (2-6)  and  (2-9)  we  then  obtain  the  potential  equation  as 


V-4  + 2k 


r>" 


(2-10) 


where  f is  a known  source  distribution  (p  ). 

source 


For  fair  weather  condition  (i.e.,  without  a charge  source) , the  potential 
equation  is  simply 


V2  + 2k 

OZ 


J«J>  = 0. 


(2-11) 


III.  Fair  Weather  Fields 

In  order  to  solve  equation  (2-11),  we  first  write  out  the  boundary 
conditions  as: 


*©• 

II 

< 

o 

at 

z = H, 

and 

<t>  = 0 

at 

z = 0. 

is  the 

upper 

boundary 

of 

potential  difference  between  the  top  of  the  atmosphere  and  the  earth 
surface . 

The  solution  for  the  potential  V at  any  altitude  in  the  atmo- 
sphere is  simply 

[ i _c-2kz  ) 


I. 

1_e-2kH  J 


(2-12) 
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The  electric  field  is  calculated  from  f = -^V  such  that 

E = 0, 
x 

E = 0 , and 

y 

V2k 

E = - . 

z h _ e-2kH 

Only  the  vertical  component  exists  and  it  is  in  the  downward  direction. 

V2k 

Defining  E = as  the  field  intensity  near  the  ground, 

0 1 - e-2kH 

we  obtain 

E = E e_2kH  (2-13) 

z 0 

the  same  result  obtained  by  Holzer  and  Saxon  (1952)  and  Mozer  (1971). 

The  average  observed  values  of  EQ  over  land  is  about  130  volts/m 
under  fair  weather  conditions.  Based  on  this  value.  Table  2-1  gives  the 
expected  field  intensities  and  potentials  with  respect  to  ground  at  dif- 
ferent altitudes.  The  Table  also  shows  the  difference  between  the  clas- 
sical picture  of  spherical  capacitor  and  the  more  recent  model  of  a leaky 
capacitor.  In  the  classical  picture,  the  field  would  fall  off  to  zero  at 
the  top  of  the  atmosphere  (50-70  km).  However,  the  new  model  allows 
0.1  m%  to  1 mV/m  of  field  at  70  km  during  fair  weather  which  is  in  the 
same  order  of  magnitude  as  fields  of  magnetospheric  origins.  During 
disturbed  weather  conditions,  more  pronounced  fields  will  be  observed. 

Figure  2-1  is  a comparison  between  calculated  values  and  the 
observed  ground  and  airborne  averages  during  fair  weather. 
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Table  2-1 

Electric  Fields  and  Potential  at  different  altitudes 
in  fair  weather  computed  from  (2-13)  and  (2-12) 
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IV.  Solution  for  Thunderstorm  Fields 

The  potential  equation  for  thunderstorm  conditions,  when  there  is 
additional  charge  source  in  the  air,  is 


V2  + 2k 


IF 


0 = - — 


(2-10) 


Tne  solution  of  equation  (2-10)  is  complicated  by  the  fact  that  the  charge 
distribution  is  spread  over  a finite  volume.  In  solving  differential 
equations  such  as  (2-10)  with  either  Dirichlet  or  Nuemann  boundary  con- 
dition the  Green's  function  method  is  both  elegant  and  effective.  The 
Green's  function  is  the  solution  of  the  potential  for  a point  source.  The 
principle  of  superposition  allows  us  to  calculate  the  potential  for  the 
actual  source  distribution. 

The  first  step  in  solving  (2-10)  is  to  generate  a Green's  function 
for  the  equation.  Assuming  azimuthal  symmetry,  the  Green's  function 
depends  on  (r,z;  r',z')  and  satisfies  the  equation 


+ 2k  h 


The  solution  of  (2-10)  is  then 


G = - -r-  6 (x-x ' ) 


(2-10') 


<Kr,z)  = 


T (volume) 


G(r,z;r',z')  f(r',z')  di ' 


in  infinite  space  with  no  boundary  surfaces. 

To  find  G we  consider  the  Fourier  transform  of  both  sides  of  (2-10'). 
The  delta  function  on  the  right  has  the  representation, 


if,,.  , iK„(z-z')  iKr  • (r-r ' 

6 (x-x  ) = dK  dd)  k d<  e z e r 

(2tt)  3 J z r r 


in  cylindrical  coordinates.  Therefore 


-24- 


G(r,z;r\z')  = fff  die  d<f>  k d<  g(<  ,ie  )eilcz(z-z')eiKr'  (r"r’>  (2-13) 

z r r r z 


The  Fourier  transform  is  to  be  determined.  Substituting  (2-13) 


into  (2-10')  gives 


die  dij)  K dk;  g (te  ,k  )(-ie,-ie2 
z r r r z r z 


* ti  \ iK,(z-z')  iKr*(r-r') 
+ 2ik  k )e  z e r 


1 ....  , iie,(z-z  ) iKr*(r-r  ) 

— ,-s , die  d<}>  k die  e z e r 

e0  (2tt)  3 z y r r 


(2-14) 


By  comparing  both  sides  of  (2-14)  we  get 


g(<r,Kz) 


e0  (271) 3 1 r z 


[ie2  + k2  - 2ik  k ]' 


(2-15) 


Substitute  (2-15)  back  into  (2-13) : 


1 fff  iKz (z-z ')i<r*(r-r') 

G(r,z;r',z')  = — .--cy  dr.  d(|)  die  ie 

C°(2lT)  iJJ  2 r r ie2  + k2  - 2ik  K 


(2-16) 


by  using  the  relation 


ei?r*(r-r')  _ eiKr|r-r'|  cos 


l iV1"  J (<|r-r'|) 


(2-17) 


which  can  be  proved  separately  (J  ()  is  the  Bessel  function  of  first  kind 

m 


of  iu^  order). 


G then  becomes 


iie,(z-z' ) 


G(r,z;r  ,z  ) = — >y  die  ie  die  , 2 . ?"  r 1 < r-r  ) (2-18) 

’ ’ ’ e.(2n)z  ! z r r (K2+Kz-2ik  k ) 0 r'  1 

o > Jr,  z r z 
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since  after  integration  over  4>  only  m = 0 term  survives. 
The  poles  of  equation  (2-18)  are  found  at 


K 

z 


i k ± /k2  + K2 
r_ 


(2-19) 


as  shown  in  Figure  2-2.  There  are  two  possible  contours,  one  in  the 
upper  half  plane  and  the  other  in  the  lower  half  plane. 

In  order  to  make  G vanish  at  z we  choose  the  contour  closing 

in  the  upper  half  plane. 


G(r,z;r ’ ,z') 
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/ K2 

14  w 


-ik{i-A+A 


2ni 

4Tr2e 

o 


k J (k  |^4,|)e-[k+/kJ+Kr](z-?,) 
r o r 1 1 


2i/k2+<2 

r 


dK 

r 


1 

4ire 

o 


foo 


0 


K J 
r 0 


(k  r-r' 
r 


. -[k+/k2+<2 ] (z-z  ) 

)e  r dK 

r 


/k2+  k2 
r 


(2-20) 


Gradshtezn  and  Ryzhik  (1965,  6.616)  gives 


J (r/m*  ^ -1)  exp(-zm')  dm'  = 


exp  (-/r^+z2) 


(2-21) 


/r2+z2 


With  proper  substitution,  (2-20)  becomes 


G(r,z;r' ,z') 


1 -k(z-z') 

47ie0 


-k/j  r-r  | 2 +(z-z'  )2 
e 

/ r-r^+Cz-z')2 


(2-22) 


Substitute  (2-22)  back  into  (2-12): 


<t>(r,z)  = 


4tte 


f(r'.z')  e 


-k( 


0 J'l  i 


-k/\r-r'\2Hz-z’f 

s-z  ) e 1 

/\r-1'f+(z -t'f 


dt ' 


(2-23) 


where  f(r',z')  is  the  source  function  and 

dx'  = r'  dr'  di}>'  dz' 

is  the  volume  element  in  cylindrical  coordinates. 

With  (2-23),  we  are  now  ready  to  solve  different  potential  problems  in  the 
atmosphere . 

The  mathematics  references  of  this  section  are:  for  Green's  function 
method,  Jackson  (1966),  Chapters  3 and  6;  for  properties  of  Bessel  functions 
and  contour  integration,  Hilderbrand  (1955),  Chapters  4 and  9 respectively. 


V . Different  Charge  Sources 

It  is  currently  accepted  that  a thundercloud  consists  of  two  portions 
of  approximately  equal  and  opposite  charges  of  electricity,  the  upper  part 
being  positively  charged  and  the  lower  part  negatively  charged.  The  charge 
separation  mechanism  (or  mechanisms)  is  not  well  understood  yet.  This 
dipolar  structure  will  change  the  electric  field  intensity  both  above  and 
below  the  cloud.  The  sign  and  magnitude  of  the  charge  depend  on  the  relative 


-27- 


position  of  the  observer  (Chalmers,  1967).  Furthermore,  there  have  also 
been  observations  and  studies  on  the  shielding  layer  around  the  cloud 
boundary  which  will  be  discussed  in  Chapter  4. 

The  simplest  approximations  for  the  charge  distribution  inside 
thunderclouds  include  (1)  point  approximation  and  (2)  thin  disc  approxi- 
mation illustrated  in  Figure  2-3.  Their  respective  range  of  validity 
will  also  be  discussed.  For  charge  distribution  of  a thick  disc,  cylinder 
or  other  shape,  the  solution  can  be  obtained  by  integration  (summation) 
over  contributions  of  points  or  thin  discs. 

The  first  step  in  solving  the  potential  and  electric  field  problem 


is  to  write  the  source  function.  For  a point  charge  at  r = 0 and 


z = z^,  the  source  function  can  be  written  as 


fi(r’’Z,)  = 2^  6(r,)  6(z’-zi) 


(2-24) 


such  that  the  potential  above  the  source  is 

2ir 


oo  no 


Vr’z)  = 4^; 


. 

qi 


2Ttr' 


z'  = 0 ^ = 0 4>’=0 

S( 


» x r / t -k(z-z')  -k/(z-z* ) 2+| r-r' | 2 

r ) o(z  -z.)  e e 

1 r' dr'd<t>'  dz’ 


J I r-r 


,|2+(z-z')2 


-k(z-z,-)  -k/r2+(z-z.- ) 2 

e 1 e l 


4ire, 


(2-25) 


/r‘!  + (z-Zi) 2 

The  boundary  condition  requires  that  4>  = 0 at  z = 0,  so  the  image 


1 


charge  is  introduced  which  gives 
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^(r.z) 


qi  -kiz-zi) 
4Tre0 


e~k/r2+(z-zi) 2 e-k^r2+(z+z^)2 

/r2+(z-z^) 2 /r2+(z+z^)2 


(2-26) 


such  that  the  vertical  and  horizontal  fields  are 


i 2 

Az 

r 

e = — - — y — 

iz  4ire.  L R 

0 n=l  n 

+ k(R  +Az  ) 
R v n n' 

n J 

exp 

-kfR  +Az  ) 
11  n n2 

(2-27) 


2 q * r 

y — 2-- 

Z „ 2 


iH  4 ire  R 

o n=l  n 


1 J 

— + k|  exp 
n ^ 


-kfR  +Az  ) 
*■  n n7 


(2-28) 


where  n = 1 and  2 refer  to  the  charge  and  image  charge  respectively 
such  that 


= q2  = qt 


Rj  = /r2+(z-zi)2. 


R2  = /r2+(z+zi)2 , 


AZj  = z-z 


Az2  = z+z^. 


The  potential  problem  for  a point  charge  source  can  also  be  solved 
by  finding  the  solution  of  Laplace's  equation  in  cylindrical  coordinates 
with  appropriate  boundary  conditions.  However,  this  method,  because  it 
does  not  apply  to  other  different  geometrical  shapes  of  charge  sources,  is 
not  as  handy  as  the  Green's  function  method.  The  procedure  is  shown  in 
the  Mathematical  Appendix  for  comparison. 

For  a dipolar  structure  of  a thundercloud  with  qj  at  z}  and  -q2 
at  z?,  the  total  field  change  is  simply  the  vector  sum  of  changes  due 
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to  individual  charges. 

<J>  = J 4^  , and 
i 1 

E - l E . 

i 

Consider  a case  of  typical  thundercloud  with  +20  Coul  at  6 km  and 
-20  Coul  at  4 km.  Table  2-2  shows  a comparison  of  the  electric  fields 
that  would  be  observed  using  this  model  with  the  property  of  variant  con- 
ductivity taken  into  consideration  and  the  "classical"  method  assuming 
uniform  conductivity  (Coulomb  field). 

It  is  shown  that  the  electric  field  above  the  source  drops  much 
faster  with  distance  as  compared  with  the  classical  method.  They  differ 
by  a factor  of  10  when  the  observer  is  24  km  overhead  above  the  cloud, 
over  100  when  39  km  above  and  1000  when  53  km  away. 

Figures  2-4  and  2-5  are  the  calculated  vertical  profiles  for  E and 

E , overhead  and  at  horizontal  distances  of  20  km,  50  km  and  100  km.  It 
H 

is  seen  that  the  electric  field  decreases  much  slower  at  larger  horizontal 
distances.  For  example,  the  vertical  fields  at  r = 0 and  r = 100  km 
differ  by  eight  orders  of  magnitude  when  observed  at  an  altitude  of  10  km 
(4  km  vertically  from  the  cloud  in  this  case),  but  at  70  km,  they  only 
differ  by  three  orders  of  magnitude.  It  is  also  shown  that  even  for  a 
moderately  charged  thundercloud  as  considered  in  this  case  (recent  works 
have  observed  up  to  several  hundred  coulombs  of  charge  inside  one  cell)  the 
change  in  the  E field  at  the  top  of  the  atmosphere  can  still  be  detected. 
It  is  in  the  order  of  millivolts  per  meter  within  a horizontal  range  of 
50  km,  comparable  to  fields  of  ionospheric  origins. 


Table  2-2 


Comparison  of  between  using  the  present  model 
((2-27)  with  (2k)  i 5km)  and  the  Coulomb  field 
calculation  ((2k)  *=&o) 
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The  same  technique  can  also  be  employed  to  find  the  potential  and 


electric  field  at  different  altitudes  for  the  case  of  thin  disc  charge 
distribution  instead  of  point  charge  source.  This  disc  distribution  of 
charge  is  similar  to  that  suggested  by  MacGorman  and  Few  (1975) . Assume 


the  r = 0 axis  is  defined  to  be  along  the  center  axis  of  the  uniformly 
charged  disc.  the  source  function  can  he  written  as 


fi(r',z')  = a'  6(z’-zi) 


(2-28) 


where  ot'  is  the  surface  charge  density  and  z^  is  the  height  of  disc. 
The  potential  can  then  be  written  as 


<t>(r,z)  = 


^=0  d)'  = 0 z'  = 0 


, x/  , . . , -k/r2+r  ' 2-2rr  'cosi(i,  + (z-zl ): 

a*6(z'-zi ) -k(z-z')  e 

4ff.0  e 

/r2+r ' 2-2rr ' cosij)  ' + (z-z ' ) 2 


r'  dr'  d<J)'  dz' 


r*  =0  <k'=o 


^ e~k(z-z'  )^-k/r2+r  '2-2rr  ,cos(|i,+(z-Zj)  : 
47rc0/r2+r ' 2-2rr  ,cos<}),  + (z-zi)2 


r’ dr’ d<J>’  (2-29) 


plus  the  image  term,  where  a is  the  radius  of  the  disc  such  that  q = va  a. 

In  order  to  compare  field  changes  due  to  point  and  disc  charge  sources, 
it  is  only  necessary  to  compare  the  vertical  fields  along  the  r = 0 axis 
where  the  difference  is  the  greatest.  For  a point  source,  from  (2-26): 


V2) 


M1  -k(z- 
4Tre0 


. -k(z-Zf)  -k(z+zj  ) 
z-z^)  e , e ^ 


(2-30) 


such  that 


E. 

1Z 


r=0 


_<L_ 

4ireo 


e-2k(z-Zl) 

(z-zi)2 


(l  + 2k(z-z^)) 


+ image  term. 


(2-31] 


For  a disc  charge  source,  from  (2-29) 


<I>i(z) 


r a 


r=0 


q^ir  g-k(z-zi) 
47re0  *Tra2 


3-k/r ’ 2+(z-z^) 2 


r*  =o 


/r ' 5+(z-zi) ■ 


r'  dr' 


+ image  term, 


2irc0a •' 


-2k(z-z^)  _ -k/a2+(z-z^) 2+(z-z^) 


. + image  term. 


(2-32) 


such  that 


E. 

q 



2e-2k(z-Zi)  _ 

( z-z . 

1 + 

-k/a2+(z-z-t ) 2+(z-z,- ) 
e 1 

1 z 

27TE0a2 

r=0 

_ 

/a2+(z-z^) 2> 

+ image  term. 


(2-33) 


Table  2-3  shows  field  changes  due  to  disc  charge  sources  as  compared 
to  point  charge  sources.  The  heights  of  different  charges  are  the  same  as 
before;  20  Coul  at  6 km  and  -20  Coul  at  4 km.  Different  radii  of  r=2  km 
and  r=5  km  for  the  disc  are  used  as  representative  sizes  of  a small  and 
medium-large  cell  respectively. 


From  these  calculations,  it  is  indicated  that  the  disc  model  is  only 
necessary  when  calculating  fields  near  the  cloud.  For  a 5 km  radius  cell, 


I 


j 
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the  difference  drops  to  10%  at  22  km  (16  km  away)  and  for  a small  cell  it 
only  differs  by  1-2%. 
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Table  2-3 

Comparison  of  field  changes  due  to  disc  charge  sources 
(2-33)  and  point  charge  sources  (2-31) 

(The  charge  distribution  is  20  Coul  at  6 km  and 
-20  Coul  at  4 km.) 
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Figure  2-1 


Computed  and  observed  dependences 
of  electric  field  strength  on  height 
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Figure  2-5 

Calculated  vertical  profile  for  horizontal  electric 
fields  at  different  horizontal  distances.  The  assumed 


source  is  a dipole  thundercloud  of  20  C at  6 km  and  - 20  C 


-42- 


| J 

r i 


CHAPTER  3 


ATMOSPHERIC  FIELDS  MAPPED  INTO  THE  IONOSPHERE 


I . Theory  and  Solution  without  Boundary  Conditions 

In  Chapter  2,  the  mapping  problem  of  atmospheric  electric  fields 
within  the  atmosphere  was  extensively  discussed.  The  upper  boundary  of 
the  atmosphere  was  assumed  to  be  at  the  altitude  of  70  km.  Yet  beyond 
that,  there  is  still  another  exciting  game  for  "atmospheric  electricians" 
to  play  instead  of  working  with  a classical  non-permissive  shielding  layer 
which  prevents  electrical  interactions  between  the  atmosphere  and  ionosphere. 

Beyond  70  km,  the  conductivity  of  air  becomes  a tensor  as  shown  in 
(1-4)  instead  of  a scalar  as  is  the  case  for  the  atmosphere.  All  the  other 
electric  equations  stay  practically  the  same.  The  basic  equations  between 
70  km  and  150  km  are  thus: 


a = 


-o. 


0 

0 


(3-1) 


(z-axis  is  along  B) 
^ x E = 0 or  E = -^<t> 


-*  “*■  r* 

ja  = o • E or  l 

i=x,y,z 


^ °ijEj 
^=x,y,z 


(3-2) 

(3-3) 


$ • 


j. 


$ • 


^ • E = — or  V2$  = - 


(3-4) 

(3-5) 
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Following  the  same  steps  as  in  Chapter  2 we  can  solve  the  potential 
problem  in  the  lower  ionosphere. 

Combination  of  (3-1),  (3-2),  (3-3)  and  (3-4)  gives 


„ fa2<t»  . 32*  . ao  32(J>  . 1 3a o 3*1  _ * t 

1 3x5  3y2  a2  3z  Oj  3z  3z 


(3-6) 


It  is  shown  from  (3-6)  that  the  field  is  independent  of  o2,  the 
Hall  conductivity,  because  of  the  antisymmetric  property  of  the  conduc- 
tivity tensor.  For  simplification,  we  can  express  oD  and  Oj  as 


a0z  = °00  6XP  (2k0Z) 


°lz  = a20  exp  (2kiz) 


(3-7) 


(3-8) 


where  the  second  subscripts  0 and  z refer  to  heights  at  and  above  the 
upper  boundary  of  the  atmosphere,  or  equivalently,  the  lower  boundary  of 
the  ionosphere.  (2kQ)-1  and  (2kj)-1  are  scale  heights  of  the  specific 
and  Pederson  conductivities  respectively.  For  better  approximation,  they 
can  be  divided  into  two  different  scale-height  regions,  one  from  70-110  km 
and  another  110-150  km. 

(3-6)  becomes 

a L-2,  + ^ + Zl  -d  +^o_«2k  — $ = ^ . t (3-9) 

i (3xz  3y2  Oj  3z^  0j  0 3zJ  3g 

Without  repeating  similar  steps  as  in  the  previous  chapter,  we  can  write 

P = • j 

rcmii*0o  fT  J e 


Pair  Go 


such  that 


_3  , 3 , On  _3  , o0  9 I » 

[to?  + w ^ 3^  + 57  2k0  aij* 


(3-10) 


and  Green's  function  satisfies 


_lL  + _9L  + £<i_iL  + ^2k  _L|g  = 

3x2  9y2  Oj  9z2  Oj  0 9zj 


6(x-x')  6(y-y')  6(z-z')  (3-11) 


The  Fourier  transformer  g should  be  such  that 


G(r,z;r\z')  = f dK  k dK  d*  g«  ,K  ) e1^2"2’*  eilcr*(^’) 


z r r 


(3-12) 


Substitute  (3-12)  into  (3-11)  and  compare  both  sides: 


dK  K d<  d<f>  g (k  *k  ) -K2  - ^ 

z r r r BV  r z'  1_  r Oj 


^ k2+  i^  2k, 


Z Oj  0 


KJ  eiKz(2-2,)eiKr*(r-r') 


1 , . - ,,  iK,(z-z  ) iK,.‘(r-r  ) 

- dK  K d<  dd)  e z e r 

e„(2tt)3  z r r y 


g(Kr’Kz)  e„  (2tt) 


HF  B + Kz  - 21  ^ ko  %]“  <3-13) 


Substitute  (3-13)  into  (3-12)  and  use  the  relation  (2-17) : 


f°°  r°°  -K_(z-z')  „ - 

! e 2 J0(Krk-r'|) 

G(r.z:r,,z')  = — . - v dK  k dK  

E.<2")*  2 ' < 

0 °i  z r °i  0 2 
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The  poles  are  at 


K 
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i(k  ± /k.?  + — K2 

l 0 ° °0 


We  choose  the  contour  closing  in  the  upper  half  plane.  Contour  integration 
then  gives 
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The  potential  is  thus 
<t>(r 


,z)  = [ G(r,z;r',z')  f(r',z')  dx  ’ 

h' 


(3-14) 


where 
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and  f is  the  source  function. 

Consider  a simple  case  of  a point  charge  source  q in  the  lower  iono- 


sphere at  r = 0 and  z = z above  70  km.  The  potential  at  altitude  z 


0 


above  the  source  is  simply 


*(r,z)  = j*-  e-k0(z-z0>  -e-°  -q' 

47T£0  


■~k  ^ -f*  (z“Z  ) ^ 

n ~ 0__  + image  term  (3-15) 


J r2  + ( z— 2 


Oj  - (z"zo>‘ 


without  any  other  boundary  conditions. 

Along  the  r = 0 axis,  (3-15)  reduces  to 


<l>  (r  = 0,  z) 


3J 

4Tie0 


e 


- k0(z~z0' 

"(z-z0) 


(3-16) 


the  same  form  as  (2-30) . 

At  all  other  distances  (r  = 0),  however,  the  potential  (and  the  field) 

decreases  much  faster  than  it  would  in  the  atmosphere  because  ^-2-  >1  in 

°i 

the  ionosphere.  and  Oj  are  the  same  at  h j<  70  km.  Above  that  O0 

increases  faster  than  O. ; at  150  km,  they  differ  by  4-5  orders  of  magni- 
tude (Figure  1-4).  Above  150  km,  the  specific  conductivity  becomes  so 
large  that  magnetic  field  lines  can  be  considered  as  equipotentials . The 
mapping  technique  developed  in  this  thesis  will  thus  be  confined  to  the 
region  below  150  km. 


II . Boundary  Conditions 


: 


Now  we  can  proceed  to  the  mapping  problem  of  atmospheric  electric  fields 
into  the  ionosphere.  We  will  assume  the  simple  point  charge  sources  because 
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it  was  shown  in  Chapter  2 that  all  the  weather-induced  charge  sources  can 
be  approximated  as  point  sources  at  large  distances. 

The  problem  yet  to  be  solved  of  field  mapping  between  the  atmosphere 
and  the  ionosphere  is  just  a boundary  value  problem  of  different  conduc- 
tivities. This  is  similar  to  boundary  value  problems  of  different  dielec- 
trics because  the  dielectric  and  conductivity  represent  the  real  and 
imaginary  parts  of  the  complex  dielectric  coefficient  respectively.  In 
order  to  complete  the  mapping  solution,  we  can  take  similar  steps  as  in 
solving  boundary  value  problems  with  different  dielectrics  which  are 
shown  in  many  electrodynamics  books  (cf.  Smythe,  1950). 

Assume  a point  charge  q embedded  at  P(z  = zQ)  in  Region  1 (atmo- 
sphere). The  boundary  surface  is  at  h(z  = 70  km)  which  separates  the 
first  medium  from  Region  2 (ionosphere)  as  shown  in  Figure  3-1.  It  is 
natural  to  locate  an  "image"  q'  in  Region  2 at  the  symmetrical  position 
P'.  The  potential  at  height  z in  Region  1 described  by  cylindrical 
coordinates  is  then 


1 Aire. 


;-k0 1 (Z+Ri ) 


rkoi(z’+RJ) 

R{ 


(3-17) 


In  Region  2,  since  there  is  no  charge  source,  it  must  be  a solution  of 
Laplace's  equation  without  singularities  in  that  region.  Clearly  the 
simplest  assumption  is  that  for  z > 70  km,  the  potential  is  equivalent 
to  that  of  a charge  q"  at  P of  the  actual  charge  st^ch  that 

- q"  e~k02(Z+R2) 

2 " 4ire0  R2 


(3-18) 
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where  (2kQj)  1 and  (2k Q2)  1 are  scale  heights  of  specific  conductivities 
in  the  atmosphere  and  ionosphere  respectively  and 

z = z - z0 

Z'  = 2h  - z„  - z 


R>  ^ + Z? 


R2=/^r2+Z2>  and 

12 


R’  = /^-L  r2  + Z'2 

1 / 0,  . 


The  boundary  conditions  required  at  z = 70  km  in  this  case  are: 
(1)  tangential  condition 


<Di  _ = *2 

lim  h ->-  70 


(3-19) 


lim  h -*■  70 


(2)  normal  condition 


a 


3£i_ 

o l 9z 


= a 


lim  h ->  70 


30? 

02  ~dT 


(3-20) 


lim  h -»■  70 


where  subscripts  ] and  2 refer  to  the  atmosphere  and  ionosphere  up  to 
150  km  respectively.  Since  the  conductivity  is  continuous  at  the  boundary, 
(3-20)  is  equivalent  to 


3<t>, 

3z 


lim  h -*•  70 


9z 


lim  h -*  70 


(3-21) 


Also  at  h (z  = 70  km),  = 1 and  Z = Z’  = 70-zQ,  R = R’  = /r2  + Z2. 
Boundary  condition  (3-19)  implies  that 

(q+q’)  e“k0l(Z+R)  = q"  e~k02(Z+R) 


(3-22) 


■ 
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Boundary  condition  (3-21)  is  reduced  to 


q-a.  = q»rc-(k,?-kn,)(z+Rrirz/R  + k02  (Z+R) 
q q q L J Lz/R  + k01(Z+R)_ 


(3-23) 


Solving  (3-22)  and  (3-23): 


q' 


= q 


(k01-k02)(Z+R) 

(k0  j+k02) (Z+R)+2Z/R 


(3-24) 


and 


q 


It 


2q 


Z/R  + knl (Z+R) 
<k01+k02)(Z+R)+2Z/R 


e t (k0 j -k0 2 ) (Z+R) ] 


(3-25) 


Substituting  them  back  into  ( 3-17)  and  (3-18) , we  obtain  potential 
expressions  both  in  the  atmosphere  and  the  lower  ionosphere  and  therefore 
can  calculate  the  vector  electric  fields. 

We  can  write  similar  expressions  if  the  lower  ionosphere  is  divided 
into  two  different  scale  height  regions.  Figure  3-2  shows  the  conductivity 
profiles  below  150  km  for  both  day  and  night.  The  scale  heights  are 
(2k02)  1 = 3 km,  (2k03)  1 = (2k12)  1 = 6 km  and  (2k13)  1 = -24  km, 

where  the  first  subscripts  0 and  1 refer  to  specific  and  Pederson 
conductives,  the  second  subsecripts  2 and  3 refer  to  regions  70-110  km 
respectively  for  nighttime  and  70-120  km  and  120-150  km  for  daytime.  Since 
the  electric  fields  are  only  affected  by  the  relative  scale  heights  of 
conductivities  but  not  their  absolute  values,  the  field  mapping  equations 
are  the  same  for  both  day  and  night. 
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III.  Results 


Electric  fields  mapped  into  the  ionosphere  can  now  be  calculated  and 
plotted  for  higher  latitude  regions  from  the  results  (3-17) , (3-18)  (3-24) 
and  (3-25).  Figure  3-3  shows  variations  in  the  vertical  fields  at  dif- 
ferent horizontal  distances.  The  vertical  field  decreases  rapidly  with 
altitude.  In  other  words,  vertical  fields  induced  by  weather  in  the 
atmosphere  dissipates  fast  in  the  ionosphere.  This  is  due  to  the  fact 
that  specific  conductivity  increases  rapidly  with  altitude.  At  even 
higher  altitudes  (from  the  F2  region  around  200  km  and  higher)  the 
conductivity  is  so  high  that  the  magnetic  field  lines  become  equipotentials 
which  force  the  electric  fields  along  them  to  go  to  zero. 

Figure  3-4  shows  variations  of  the  horizontal  fields.  Within  the 
50  km  horizontal  range  considered,  they  vary  in  the  same  bulk  order  of 
magnitude  as  the  vertical  fields.  However,  at  higher  altitudes,  200  km 
or  above,  when  o0/Oj  >>  1 (for  example,  at  200  km,  °0/oi  % 10'’-106 
and  at  500  km,  O0/O(  ^ 108)  and  such  that  /o0 /O^  • r >>  z,  the  vertical 
field  vanishes  and  the  horizontal  field  approaches  an  asymptotic  value. 

In  other  words,  E^  becomes  a function  of  r,  independent  of  z. 

To  demonstrate  the  case  mathematically,  we  start  from  (3-18) 


e-k0(Z+R) 


(3-18) 


where  z = z-zr 


R = r2+z2 

/oi 


or  Lim  R2  = 


such  that 


This  leads  to 


Figure  3-5  is  a diagram  of  different  conductivity  profiles  in  the 


ionosphere  up  to  1000  km.  It  is  shown  that  above  the  F2  region  ('v250  km) 


O . remains  constant  which  implies  that  k -*■  00 . Thus  (3-26)  and  (3-28) 


are  simply 


Lim  $ = 


o,  1 


4 Tie. 


o r 


(3-29) 


and 


Lim  E 


= sl! — /2j 

t “ rt 


r 4Tie0  °o  r 


(3-30) 


It  is  also  of  interest  to  observe  the  corresponding  horizontal 


distance  for  maximum  E at  different  altitudes.  It  is  a simple  relation- 


ship for  Coulomb  field  conditions.  Since 


(ji)  oc 


(r2+z2)1/2 


E « — rr 

r ' ,2  ,2,  Vz 


(rz+z^) 


__ 


and 


such  that 


= 0 implies 


2r2  - z2  = 0 


f=  z = 0.7z. 


However,  in  the  atmosphere-ionosphere  system,  the  situation  is 

quite  different.  In  the  isotropic  atmosphere,  the  ratio  r|  ]z 

br  max  / 

starts  from  ^0.2  at  lower  altitude  and  increases  to  Ml.  5 at  20  km, 
then  drops  slowly  till  the  top  of  the  atmosphere.  In  the  anisotropic 
ionosphere  the  ratio  drops  at  a much  faster  rate  with  increasing  altitude. 
At  120  km,  E is  at  r 'v  0.5  km  which  gives  the  ratio 

1 IllaX 

r | ^ I z 0.004.  Table  3-1  gives  in  more  detail  the  horizontal 

r max  / 

fields  at  different  altitudes.  We  have  again  considered  the  same  charge 
source  as  in  the  previous  sections:  +20  C at  6 km  and  -20  C at  4 km. 

At  lower  latitude  regions  where  magnetic  field  lines  are  no  longer 
along  z axis,  the  differential  equation  to  be  solved  is  much  more  com- 
plex. The  conductivity  tensor  is  of  a more  general  form 


f OjSin  I + cos  I 02sinl  sinl  cosl 


— o 2 sinl 


0 2 cosl 


(a.-Oj)  sinl  cosl  -a 2 cosl  0j  cos2I  + oQ  sin2I 


where  I is  the  magnetic  dip  angle;  at  pole  I = 90°  and  at  equator  1=0 
for  a dipole  model. 


L 
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At  equator,  it  is  reduced  to 


0 -a. 


such  that  the  potential  equation  to  be  solved  is 


iff  + iff  + iff  _ 2k2  ii  + 2kl  i*  _ _ _L 

Oi  8?  37  3?  K2  o,  3y  9z  e0 


(3-31) 


Difficulty  arises  because  equation  (3-31)  does  not  have  the  property  of 
spherical  symmetry.  Different  approximation  methods  and  computer  tech- 
niques will  have  to  be  used  in  order  to  obtain  the  electric  field  pro- 
file. Furthermore,  in  the  upper  atmosphere  of  lower  latitude  regions, 
tidal  winds  interact  with  the  geomagnetic  field  and  charged  particles 
which  induces  electric  fields  and  current  flow.  According  to  this 
dynamo  theory  (first  developed  by  Chapman  and  Bartels,  1940),  the  in- 
duced electric  fields  bear  both  solar  and  lunar  variations  which  will 
add  more  difficulty  on  filtering  out  the  thunderstorm  electric  fields  in 
the  lower  latitude  ionosphere.  Generally  speaking,  however,  we  can 
anticipate  that  at  equatorial  region,  the  component  of  the  electric  field 
parallel  to  the  magnetic  field  lines  will  dissipate  faster  than  other 
components,  and  E will  reach  an  asymptotic  value  at  higher  altitudes. 
Furthermore,  since  ionospheric  conductivities  are  much  higher  at  equa- 
torial regions,  especially  during  the  day,  than  at  higher  latitude 
regions,  the  electric  fields  will  map  less  efficiently  in  equatorial  than 
higher  latitude  regions.  In  short,  electric  fields  induced  by  thunder- 
storms in  the  equatorial  region  are  less  efficient  in  being  mapped  into 
the  lower  ionosphere  than  those  produced  at  higher  latitude  regions. 
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Figure  3-2a 

Conductivity  profile  approximating  the  daytime  condition 
above  70  km  and  the  Cole  and  Pierce  profile  below  70  km. 
The  numbers  next  to  the  curves  are  scale  heights  in  kilo- 
meters. (Park  and  Dejnakarintra,  197^) 
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Figure  3-2b 

Conductivity  profile  approximating  the  nighttime  condition 
(Park  and  Dejnakarintra,  197^) 


Figure  3-3 

Vortical  fields  mapped  into  the  ionosphere 


Figure  3-4 

Horizontal  fields  mapped  into  the  lower  ionosphere 
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Figure  3-4 


f I 
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Horizontal  electric  fields  mapped  into  the  lower  ionosphere. 
The  source  is  an  electric  dipole  of  20  G at  6 km  and  -20  C 
at  4 km. 
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Figure  3-5 

Large  scale  conductivity  profiles  in  the  ionosphere. 

& pi  4 <40  are  pederson,  Hall  and  specific  conductivities 

respectively. 

A refers  to  the  region  within  aurora  arc.  C refers  to  the 
region  of  the  undisturbed  ionosphere  surrounding  the  arc. 
(Bostrom,  1964) 
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CHAPTER  4 

COMPARISON  BETWEEN  THE  MODEL  AND  DATA  FROM  A BALLOON  FLIGHT 

I . General  Description  of  the  Balloon  Flight  and  Data 

In  August  1969,  a total  of  24  balloons  were  flown  in  four  different 
launches  from  six  Canadian  sites  (Figure  4-1)  with  each  of  the  balloons 
in  a group  simultaneously  returning  data  for  up  to  one  day.  Approxi- 
mately 370  hours  of  data  are  available  and  the  vector  electric  fields 
can  be  calculated.  The  sensitivities  of  the  vector  components  were  such 
that  during  fair  weather  conditions,  the  vector  electric  field  was 
recorded  at  less  than  half  of  the  full  scale  in  that  channel,  while  the 
slowly  time-varying  horizontal  fields  (often  associated  with  magneto- 
spheric  activity)  were  displayed  in  the  two  other  channels.  Thus  fair 
weather  conditions  can  be  studied  in  detail.  During  several  flights 
where  balloons  passed  over  thunderstorms,  atmospheric  fields  were  strong 
enough  at  times  to  saturate  the  continuous  channels  for  the  three  com- 
ponents. However,  large  dynamic  range  electric  fields  signals  were  also 
transmitted  in  commutated  forms. 

The  balloon  payloads  have  been  described  in  various  papers  (Mozer 
and  Serlin,  1969;  Mozer  and  Manka,  1971;  and  Mozer,  1971).  A typical 
structure  is  as  follows: 

The  balloon  was  separated  from  the  payload  by  a rope  of  about  50 
meters  long  such  that  its  presence  would  not  distort  the  fields  at  the 


1 
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payload.  Volumes  of  different  balloons  ranged  from  about  1600  cubic 
meters  to  around  7100  cubic  meters  such  that  ceiling  altitudes  of  26  km 
to  34  km  were  reached. 

The  electric  field  sensors  were  three  mutually  orthogonal  pairs  of 
dipole  probes,  two  horizontal  and  one  vertical.  Each  dipole  probe  pair 
consisted  of  an  insulated  ( 'V  1014  ohms)  rod  about  3 meters  long  with  a 
conducting  plate  at  each  end.  Each  electric  field  component  was  deter- 
mined as  the  potential  difference  between  the  appropriate  pair  of  plates 
divided  by  their  separation  distance.  The  payload  was  driven  by  an  elec- 
tric motor  such  that  it  rotated  with  a period  of  20  to  25  seconds.  An 
on-board  magnetometer  provided  orientation  information  which  allowed  reso- 
lution of  the  horizontal  field  into  components  in  the  magnetic  north  and 
east  directions. 

In  order  to  measure  the  ambient  electric  field  with  accuracy,  it 

was  necessary  to  utilize  very  high  input  impedance  preamplifier  connected 

to  the  sensors.  The  signal  of  interest,  V,  causes  a current  of 

V/R  . to  flow  through  the  electronic  input  resistance,  R . . . 

circuit  6 circuit 

This  current  returns  through  the  atmosphere  along  a path  with  resistance 

R . Unless  R . . „ » R , an  attenuated  signal  would  be  measured 

atm  circuit  atm 

because  an  appreciable  amount  of  V would  then  appear  across  Ratln*  T^e 

value  of  R depends  on  the  conductivity  of  air  at  different  altitudes, 
atm 

For  the  balloon  system  considered,  R ~ 1013  ohms  near  the  surface  of 

atm 

the  earth  and  R ^ = 1011  ohms  at  30  km.  Thus  an  input  resistance 

atm 

greater  than  1014  ohms  was  required  for  field  measurements.  Electric 
circuits  meeting  this  requirement  were  flown  on  the  balloons.  The  four 
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subcarrier  inputs  to  the  FM  - FM  telemetry  transmitter  consisted  of  the 
three  components  of  electric  fields  and  a commutator  through  which  other 
ranges  of  sensitivity  on  the  electric  field  measurements,  magnetometer 
and  horizontal  orientation  information,  housekeeping  temperature,  vol- 
tages and  currents  were  transmitted. 

In  this  thesis,  the  effort  has  been  concentrated  on  the  August  3, 
1969,  balloon  flight  launched  from  Penhold,  Canada,  at  0328  UT  (2028  LT) . 
Data  were  collected  for  about  21  hours.  They  included  information  on 
thunderstorm  activity  in  the  vicinity  during  the  first  portion  of  the 
flight  and  also  the  recovered  fair  weather  field  in  the  later  portion. 

The  following  pieces  of  data  (and  reduced  data)  have  been  obtained: 

(a)  continuous  raw  data  of  the  whole  flight  showing  saturation 
of  field  during  the  thunderstorm,  spikes  due  to  lightning 
flashes  and  the  continuous  fair  weather  field  data  after 
recovery  of  the  storm; 

(b)  plots  of  vector  fields  with  60  minute  averages  on  compressed 
scale  (0.4  in/hr)  of  the  whole  flight; 

(c)  plots  of  vector  fields  with  15  minute  averages  on  the 
same  scale  as  (b)  of  the  whole  flight; 

(d)  computer  printout  and  plots  of  vector  fields  with  30  second 
averages  on  expanded  scale  (12  in/hr)  from  0500  UT  to  1000 
UT  for  d-c  thunderstorm  electric  field  study;  and 


(e) 


computer  printout  and  plots  of  electric  fields  with 


0.5 


i 


second  averages  on  more  expanded  scale  (4  in/hr)  from 
0800  UT  to  0850  UT,  when  lightning  activity  was  present, 
for  detailed  study  of  lightning  induced  fields. 


If 
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The  thunderstorm  passed  Penhold  area  at  around  0300  UT  with  a velo- 
city of  15  miles/hr  to  the  east.  The  balloon,  launched  at  0328  UT,  was 
first  drifting  to  the  east  until  around  0630  UT.  It  then  moved  to  the 
west  at  15  to  20  miles/hr  and  was  overhead  at  Penhold  again  around  0700 
UT.  The  ceiling  height  was  around  28.7  km  (94,000  ft). 

Figure  4-2  shows  the  raw  data  format  of  balloon  field  data.  Figures 
4-3a  and  b are  portions  of  raw  data  indicating  the  presence  of  thunder- 
storm in  the  vicinity  and  distinct  lightning  events  respectively. 

II . Fair  Weather  Field  Data 

Figure  4-4  shows  the  fair  weather  vector  electric  fields  during  the 
later  portion  of  the  flight.  They  can  be  discussed  separately  as  the 
vertical  and  horizontal  components  of  the  electric  field. 

(1)  Vertical  field: 

The  average  reading  for  the  vertical  field  was  at  300  mV/m. 

This  agrees  with  the  prediction  of  (2-13)  for  z at  28  km,  the  balloon 
altitude  at  the  time,  with  scale  height  of  conductivity  (2k)-1  = 5 km. 

We  have  to  keep  in  mind  this  downward  field  as  the  background  value  when 
analyzing  field  data  during  the  thunderstorm  period. 

(2)  Horizontal  field: 

The  fair  weather  horizontal  fields  observed  by  balloon  measure- 
ments in  higher  latitude  regions  are  mostly  of  magnetospheric  origins 
(Mozer  and  Manka,  1971).  Enhancement  in  horizontal  electric  field  with 
increasing  Kp  (a  number  ranging  from  0 to  9 indicating  the  degree  of 
disturbance  of  the  geomagnetic  field)  has  largely  been  observed.  During 
quiet  conditions,  the  observations  can  be  best  explained  by  the  sunward 
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E * B convection  model.  An  overall  average  of  32  balloon  flights  in 
the  aurora  zone  indicated  that  the  dominating  field  component  is  south- 
ward near  dawn  (0600  LT)  and  northward  near  dusk  (1800  LT)  (Mozer  and 
Lucht,  1974).  They  also  showed  that  for  the  south-north  field  component 
the  maximum  field,  which  occurs  at  dawn  or  dusk,  is  less  than  25  mV/m 
and  it  gets  close  to  zero  at  0000  LT  and  1200  LT  during  periods  of  low 
Kp.  During  times  of  high  Kp,  it  fluctuates  between  ± 40  mV/m.  For 
west-east  field  component,  it  does  not  get  beyond  a few  mV/m  during  quiet 
times  and  can  go  up  to  25  mV/m  during  disturbed  times. 

At  around  0900  UT  August,  1969,  the  Kp  index  was  at  2-  (Solar- 
Geophysical  Data,  1969)  which  indicated  a relatively  quiet  condition. 

As  a result,  the  west-east  field  component  should  be  zero;  and  the  south- 
north  component,  since  it  was  near  0200  LT,  should  also  vanish.  This 
agrees  with  observations  in  Figure  4-4.  Furthermore,  we  should  also  bear 
in  mind  that  during  the  thunderstorm  period  around  0600  UT  (2300  LT), 
there  should  not  be  much  significant  field  of  magnetospheric  origins  so 
that  the  electric  field  observed  during  that  period  would  be  all  weather- 
induced  . 


atmosphere  can  be  predicted  from  calculations  of  (2-26)  and  (2-27).  During 


the  Penhold  balloon  flight  of  August  3,  1969,  observations  of  thunderstorm 
measurements  were  between  0500  UT  and  0630  UT.  At  this  time,  the  balloon 
was  moving  to  the  east  relative  to  the  storm  at  about  20  miles  per  hour. 
According  to  nearby  radiosonde  data  from  Edmington,  70  miles  northwest  of 
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Penhold,  the  freezing  level  was  at  around  lh  km  altitude  at  the  time 
(Monthly  Bulletin,  Canadian  Upper  Air  Data,  1969) . It  is  generally 
accepted  that  the  negative  charge  region  inside  a thundercloud  is  concen- 
trated near  the  -5°C  isotherm.  The  ceiling  height  of  the  clouds  was 
around  7 km.  It  is  then  reasonable  to  assume  a dipole  cloud  model  with 
the  centers  of  negative  and  positive  charge  regions  at  4 km  and  6 km  res- 
pectively. Since  the  balloon  was  at  sufficiently  high  altitude  comparing 
to  the  cloud,  point  charge  assumption  is  a good  enough  approximation  for 
different  kinds  of  models;  for  example,  cylindrical,  disc  models  either 
upright  or  tilted. 

Also  according  to  the  same  radiosonde  information,  at  balloon  alti- 
tude, the  temperature  was  at  -40°C,  pressure  at  20  mb  and  the  wind  was 
easterly  at  6 miles/hr.  And  at  the  altitude  near  the  storm,  the  tempera- 
ture was  -15°C  along  the  500  mb  isobar  and  the  wind  was  southwest-westerly 
at  around  14  miles/hr. 

By  comparison  between  the  maximum  changes  in  vertical  and  horizontal 
fields,  the  charge  concentrations  were  estimated  at  ± 20  Coul.  A Carte- 
sian coordinate  system  can  be  defined  such  that  x-axis  is  to  the  east, 
y-axis  to  the  north  and  z-axis  in  the  upward  direction.  Since  the  balloon 
was  movoing  to  the  east  relative  to  the  storm,  the  horizontal  field  com- 
ponents along  and  perpendicular  to  the  balloon  trajectory  can  be  referred 
to  as  the  x and  y components  of  the  total  field  respectively.  The 
vertical  field  is  the  z component. 

Figure  4-5  shows  the  expected  vertical  field  at  balloon  altitude 
( 'v  28  km)  along  different  horizontal  trajectories  of  y = 0,  10  and  20km. 
The  field  above  the  axis  refers  to  upward  direction.  The  maximum  verti- 
cal field  changes  for  this  case  are  5.00  V/m  when  overhead,  3.44  V/m 


along  the  y = 10  km  trajectory  and  1.37  V/m  if  y = 20  km.  From  the 
figure  it  shows  the  field  decreases  fast  with  increasing  distance  in  x. 

At  x = 40  km,  it  approaches  zero.  Actually  for  a dipole  source,  the 
vertical  field  would  change  its  sign  at  some  horizontal  distance.  Closer 
than  this  distance  the  vertical  field  is  upward  and  beyond  this  distance 
the  field  is  downward  (the  same  direction  as  the  fair  weather  field). 
However,  at  balloon  altitude,  the  field  does  not  change  its  sign  until 
x > 100  km.  Since  the  magnitude  of  the  field  change  is  already  very 
small,  this  effect  is  not  at  all  noticeable  and  would  not  be  detected 
by  the  balloon  field  sensors. 

Figure  4-6  is  the  expected  electric  field  changes  in  the  horizontal 
plane  along  the  trajectory  (E^) . For  the  case  that  the  balloon  moves 
to  the  east  relative  to  the  storm,  the  direction  of  this  component  of 
electric  field  is  such  that  when  the  balloon  is  approaching  the  thunder- 
cloud (negative  x)  the  field  is  to  the  west,  and  when  the  balloon  is 
receding  from  the  storm,  the  field  is  to  the  east.  The  maximum  hori- 
zontal field  changes  along  the  trajectory  due  to  the  thundercloud  con- 
sidered are  0.86  V/m  at  x = 10  km  along  the  y = 0 trajectory,  0.60  V/m 

at  x = 12  km  if  y = 10  km  and  0.25  V/m  at  x = 15  km  if  y = 20  km. 

Figure  4-7  is  the  expected  electric  field  changes  in  the  horizontal 
plane  perpendicular  to  the  trajectory  (E  ).  The  direction  of  this  field 
component  is  such  that  if  the  trajectory  is  north  to  the  thundercloud,  a 
northward  field  will  be  experienced  all  the  way.  For  the  case  when  the 
balloon  is  flying  along  y = 0,  no  E^  will  be  experienced.  Otherwise 

the  maximum  occurs  at  x = 0.  The  values  are  0.85  V/m  when  y = 10  km 
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Figure  4-8  shows  the  averaged  field  data  in  the  vertical,  east- 
west  and  north-south  directions  between  0500  UT  and  0630  UT.  In  order 
to  trace  the  number  of  cells  encountered,  to  estimate  the  magnitude  of 
charge  involved  and  to  locate  the  relative  position  of  the  storm  with 
respect  to  the  balloon  trajectory,  the  following  steps  are  taken: 

(1)  By  studying  the  peaks  in  the  vertical  field,  we  can  deter- 
mine the  number  of  cells  with  charge  concentrations  and 
their  distances  from  each  other.  In  Figure  4-8,  it  is 
indicated  that  there  were  two  major  cells  encountered  by 
the  balloon  at  around  0535  UT  and  0605  UT  respectively; 
another  one  right  before  0500  UT  and  one  after  0630  UT. 

Since  according  to  radar  tracking  at  the  time,  the  rela- 
tive velocity  of  the  balloon  was  about  20  miles/hr  rela- 
tive to  the  storm,  it  implies  that  the  cells  were  approxi- 
mately 10-15  miles  (16-20  km)  from  each  other.  This  agrees 
with  weather  observations  which  indicated  the  average 
separation  between  cells  ranged  from  10-20  miles.  Further- 
more, the  difference  between  the  peak  vertical  field 
value  and  the  lowest  value  should  be  the  maximum  change 

of  the  vertical  field.  We  thus  obtain  the  maximum 
vertical  field  due  to  thundercloud  structure  to  be  around 
1.2  V/m  in  the  upward  direction. 

(2)  Since  the  balloon  trajectory  relative  to  the  storm  was 
along  the  east-west  direction,  the  north-south  field, 
the  horizontal  field  perpendicular  to  the  trajectory, 
should  vary  in  the  same  phase  as  the  vertical  field.  As 
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shown  in  Figure  4-8,  the  peak  readings  in  the  southward 
direction  occurred  at  generally  the  same  time  as  the 
vertical  field  which  indicated  that  the  balloon  trajec- 
tory was  south  to  the  thunderclouds.  The  only  exception 
is  the  peak  of  the  southward  field  at  0510  UT  which  indi- 
cated that  either  the  cell  that  the  balloon  just  flew 
by  had  a tilted  dipole  structure  or  the  balloon  trajec- 
tory with  respect  to  the  cell  was  diverted  from  the 
assumed  south  direction.  The  difference  between  the 
maximum  magnitude  and  the  lowest  value  was  around  0.45 
V/m.  We  can  then  estimate,  together  with  the  vertical 
field,  how  far  away  the  balloon  trajectory  was  from  the 
storm  and  the  charge  concentration  in  each  cell.  Com- 
paring results  from  Figure  4-5  and  Figure  4-7,  we  can 
conclude  that  the  charge  concentration  was  around  ± 20 
Coul  in  each  major  cell  and  the  trajectory  was  20-22  km 
south  to  the  thunderclouds. 

(3)  One  major  problem  left  is  the  almost  constant  bias  in 
the  measured  vertical  and  north-south  fields.  After 
subtracting  the  0.3  V/m  normal  downward  field,  the  elec- 
tric field  still  has  an  additional  0.3  V/m  in  the  down- 
ward direction  and  around  0.15  V/m  in  the  northward 
direction.  They  are  about  1/4  to  1/3  of  the  maximum  of 
the  field  component  in  the  opposite  directions.  It  seems 
most  probable  that  there  must  be  a negative  charge  source 
spread  over  the  thunderclouds.  The  possible  explanation 


for  this  phenomenon  is  as  follows: 

There  exists  a "shielding  layer"  in  the  boundary  region 
of  the  cloud.  The  formation  of  this  layer  is  due  to  the 
differences  in  conductivities  between  the  air  and  the 
cloud  at  cloud  boundary.  Measurements  have  indicated 
that  the  cloud  is  a much  less  conducting  medium  than  air. 
Because  of  the  dipole  structure  of  the  thundercloud,  there 
will  be  a layer  of  negative  charge  accumulated  on  top  of 
the  cloud  to  match  the  discontinuous  boundary  condition. 
There  should  also  be  a layer  of  positive  charge  near  the 
bottom.  Phillips  (1967)  has  made  theoretical  calculations 
using  the  method  of  ionic  equilibrium  at  the  cloud  boundary. 
He  showed  that  the  loss  by  conduction  to  cloud  droplets 
due  to  the  shielding  layer  charge  distribution  was  about 
1/3  that  would  occur  in  the  absence  of  such  a layer.  The 
thickness  of  the  shielding  layer  was  calculated  to  be  in 
the  order  of  tens  of  meters.  The  total  shielding  charge 
in  this  case  is  estimated  to  be  ± 5 Coul  at  the  top  and 

bottom  boundaries  of  each  cell. 

The  field  detected  by  the  field  sensors  thus  include 
(1)  normal  background  fair  weather  field  (downward),  (2) 
the  field  due  to  thunderclouds  (upward)  and  (3)  the  shield- 
ing field  (downward).  Figure  4-9  shows  how  these  fields 
combine  to  give  the  total  vertical  fields  and  how  they 
compare  with  the  observed  values.  The  horizontal  fields 
can  be  done  similarly. 

(4)  From  (1)  and  (2)  we  estimated  the  charge  concentrations 
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to  be  ± 20  Coul  in  each  thundercloud  and  the  balloon 


trajectory  was  around  22  km  south  of  the  storm.  We 


can  then  get  the  resultant  east-west  horizontal  field. 


the  field  component  along  the  direction  of  the  trajec- 


tory, by  superposition  of  the  east-west  field  due  to 


different  cells  (Figure  4-10) . The  calculated  results 


agrees  with  the  shape  of  the  experienced  field  except 


most  of  the  time  there  existed  an  unexplained  additional 


eastward  field  of  0.005  to  0.008  V/m  in  magnitude.  This 


was  not  due  to  the  large  scale  shielding  effect  just 


discussed  because  the  component  along  the  trajectory 


should  cancel  out.  It  could  not  be  of  ionospheric 


origin  either  because  the  ionospheric  condition  was 


relatively  quiet  at  the  time  and  0.08  V/m  is  too  large 


a field  compared  to  the  normal  0.02  V/m  ionospheric 


fluctuations.  It  is  still  most  probably  of  atmospheric 


origin,  namely  a general  disturbed  area  to  the  west  as 


a continuation  of  the  storm.  This  source,  if  existed, 


had  to  be  far  off  enough  such  that  peaks  in  fields  due 


to  different  cells  would  not  be  experienced  at  the  balloon. 


In  general,  the  data  agree  with  the  predictions  and  the  assumed 


models.  We  are  able  to  locate  the  storm  position  relative  to  the 


balloon  and  estimate  the  charge  distribution.  The  major  difficulties 


in  data  analyzing  include: 


(1)  Lack  of  continuous  storm  informatic  ' from  radar.  For  example. 


information  on  the  time  history  of  che  storms  and  locations  of  different 
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thunderclouds  at  all  times  were  not  available. 

(2)  Saturation  of  data  information.  In  Figure  4-8,  part  of  the 
vertical  and  north-south  field  readings  were  missing  because  of  instru- 
ment saturation.  Some  of  the  peak  values  could  only  be  approximated 
from  the  slopes  of  the  general  envelop.  As  a result,  only  estimates 
could  be  made  instead  of  exact  calculations. 

Nonetheless,  this  balloon  flight  has  been  one  of  the  very  few 
vector  electric  field  measurements  above  thunderstorms  that  indicated 
definite  field  variations  and  provided  valuable  information  on  thunder- 
cloud mapping  into  the  upper  atmosphere.  The  improvements  and  potential 
of  field  measurements  by  balloons  will  be  discussed  in  the  next  chapter. 


r 
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I V . A Study  of  the  Lightning  Events 

In  the  preceding  section  we  have  examined  the  thunderstorm  induced 
electric  fields  due  to  static  or  slow  varying  distributions  of  charges 
inside  thunderclouds.  During  a later  portion  of  the  flight  (around  0830 
UT)  the  balloon  experienced  some  lightning  activity.  In  this  section 
we  will  study  these  lightning  induced  fields  in  some  detail. 

In  a lightning  discharge,  there  is  a considerable  variation  in 

1 

charge  and  current  in  a relatively  short  time  interval.  Thus  in 
addition  to  the  electrostatic  field,  we  also  have  to  look  into  the  time- 
varying  electromagnetic  fields  associated  with  lightning  discharge. 

For  a system  of  source  of  charges  and  fields  varying  with  time, 
we  can  assume  the  time  variation  factor  to  be  e-*wt  such  that 

-iwt 

p(x,t)  = p(x)e 

and  j(x,t)  = j (x)e  iwt 


( 
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where  the  real  part  of  the  expression  is  to  be  taken  to  obtain  physical 
quantities.  The  electromagnetic  potentials  and  fields  are  also  assumed 
to  have  the  same  time  dependence.  In  the  limit  that  source  region  is 
of  a much  smaller  scale  compared  to  one  wavelength  X = 27rc/u>,  spatial 
regions  of  interest  can  be  divided  into  three  zones: 
the  near  (static)  zone:  R « X, 

the  intermediate  (induction)  zone:  R ^ X, 

the  far  (radiation)  zone:  R » X. 

In  the  near  zone  the  fields  have  the  character  of  static  fields.  In 
the  far  zone,  the  fields  are  transverse  to  the  radius  vector  and  fall 
off  as  R”1,  typical  of  radiation  fields. 

To  solve  the  time-dependent  problem,  it  is  no  longer  valid  to 
find  an  equation  for  the  scalar  potential  4>,  defined  as  E = - V4> 
because  ^ x E = - 9B/9t  = iwB  f 0.  We  have  to  define  a vector  poten- 
tial  A such  that 


and 


B = $ x t 

C dt 


The  wave  equation  for  A in  our  case  is 


(V2  + 2k  -)  A 


1 92A 
c7  9t2 


The  magnetic  induction  is  given  by 

B = $ x A 

and,  outside  the  source,  the  electric  field  is 

E = — ? x B ■ 

0) 


(4-1) 
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For  a simpler  case  of  Coulomb  field  ((2k)-1  -*•  °°)  (4-1)  reduced  to  the 

familiar  form 


V A - 


1 92A 


c2  9t2 


ce, 


For  a-  electric  dipole  with  the  dipole  moment  of 

x'  p(x')  d3x'. 


P = 


The  solution  of  the  fields  are  (Jackson,  1952) 


B = 


0) 


(n  x P)  V - (i  - -±-) 


0) 


R '*  i k R' 
i§R 


e = ~r  (n  x p)  x n 


R 


, 1 i . 1cR 


+ [3n(rrp)  - p](pT  - ^y)« 


(4-2) 


In  the  near  zone,  R « X = 2ttc/o), 


E = 3n[(n*p)  - p] 


the  electric  field  is  just  the  static  electric  dipole  field.  In  the 
far  zone,  R » X, 


2 rrR 

-+  uv  -*■  e c 

E = -^2  (n  x p)  x n — £— 


the  electric  field  shows  the  typical  behavior  of  radiation  fields. 

We  can  now  compare  the  magnitudes  of  radiation  and  static  fields 
in  the  far  zone.  From  (4-2),  we  obtain 


"radiation 


"static 


J2  2 2TT2f2R2 


2c/ 
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At  a distance  of  100  km  away  from  the  source,  the  ratio  is  about 
2 x 10~6  x f2  where  f is  the  frequency.  In  other  words,  if  the 
frequency  is  close  to  1 kHz,  radiation  and  static  fields  are  in  the  same 
order  of  magnitude  at  this  distance.  Above  this  frequency,  radiation 
field  will  dominate. 

The  electrostatic  field  is  the  field  that  we  obtained  in  Chapters 
2 and  3;  it  is  the  only  non-zero  field  before  and  after  a lightning 
charge.  The  radiation  field  represents  the  energy  propagated  away  from 
the  source  at  the  speed  of  light. 

In  the  atmosphere,  within  a range  of  around  30  km  from  the  lightning 
source,  the  dominating  field  is  electrostatic.  Beyond  this  distance, 
the  induction  and  radiation  terms  start  to  contribute  to  the  total 
field.  In  the  far  zone,  the  radiation  term  will  be  the  only  significant 
field. 


lightning  region  relative  to  the  balloon  and  to  estimate  the  average 
discharge  involved  in  each  flash.  We  first  summarize  the  data  as 
follows : 

(1)  The  maximum  change  in  vertical  field  caused  by  lightning 
was  ^0.60  - 0.70  V/m  in  the  downward  direction. 

(2)  The  maximum  horizontal  fields  were  % 0.20  V/m  and  0.25 
V/m  in  the  east  and  north  directions  respectively. 

(3)  The  ratio  between  vertical,  east  and  north  fields  was 
about  the  same  for  all  flashes.  Furthermore,  the  values 
of  fields  were  also  in  the  same  order  of  magnitude  from 
one  flash  to  another. 

The  direction  of  change  in  the  vertical  field  indicated  the  light- 
ning flashes  experienced  by  the  balloon  could  not  be  cloud-to-ground 
or  cloud-to-sky ; they  had  to  be  intracloud  which  involve  downward 
positive  discharge.  The  directions  of  changes  in  vertical  fields  due 
to  various  types  of  flashes  are  shown  in  Figure  4-12. 

Comparison  of  (1)  and  (2)  with  results  from  (2-26)  and  (2-27) 
indicates  that  the  source  region  was  still  north  of  the  balloon  tra- 
jectory as  was  the  case  earlier.  The  average  amount  of  discharge  was 
around  20  Coul  and  the  flashes  were  about  18  km  east  ahead  of  the 
balloon;  making  the  lightning  region  approximately  27  km  northeast  of 
the  balloon.  This  agrees  with  the  local  weather  observation  at  the 
time.  From  (3)  we  know  the  lightning  fields  measured  for  the  balloon 
all  came  from  the  same  general  area.  In  other  words,  there  was  only 
this  particular  lightning  active  region  within  the  range  of  balloon 
detection. 


There  were  seven  distinct  lightning  flashes  observed  by  the  balloon 

field  sensors.  The  detail  in  the  vertical  field  during  each  flash  is 

shown  in  Figure  4-13.  Table  4-1  is  a summary  of  the  amount  of  discharge, 

the  time  constant  for  field  recovery  (time  T required  for  the  field  to 

go  back  to  1/e  of  the  peak  value  Ae  ) . 

max 

From  Figure  4-13  and  Table  4-1,  the  amount  of  discharge  for  indivi- 
dual flash  is  between  18.5  Coul  and  24  Coul.  The  time  constant  for  field 
recovery  varies  between  2.5  and  7.5  seconds. 

The  problem  of  electric  field  recovery  after  lightning  has  been 
studied  to  some  extent  (Illingworth,  1971,  1972;  Vonnegut  et  al.,  1966). 
Experimental  measurements  and  statistical  analysis  of  recovery  curves 
have  been  of  great  diversity  and  non-conclusive.  However,  the  basics 
of  physics  of  field  response  to  the  charge  regeneration  at  different 
altitudes  within  the  electrostatic  range  in  the  atmosphere  can  be  out- 
lined easily.  Equations  (2-1),  (2-2)  and  (2-3)  still  hold. 


C?z  = 0O  exp  (2kz) 


x E = 0 or  E 


j = 0E 
Ja 


= - ^ 


(2-1) 


(2-2) 


(2-3) 


(2-4)  is  no  longer  valid  since  non-steady  state  conditions  have  to  be 


taken  into  consideration. 


such  that 


V (j  + e„  H ) - 0 


Table  4-1 

A general  feature  of  the  different  lightning  flashes 
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where 


j = j + j . 
a s 


Equation  (2-4)  can  tlicn  be-  rewritten  as 


3E 


^ • i + £ = — V • 1 

Ja  eo  3t  V Js 


(4-3) 


Equations  (2-1),  (2-2),  (2-3),  (4-3)  together  with  Poisson's  equation 

V2$  = - -£-  (2-5) 

eo 

then  give 


V2<I>(x,t)  + 2k  3j>^X-’t-)-  + A (V2<Kx,t))  = - f (x,t) 

o Z O o t 6 » 


(4-4) 


where  f is  the  source  function. 


At  t = 0,  V2$  can  be  derived  from  (2-26),  3/3t(72$)  can  be 


predicted  as  of  the  form 


; 


A (v2$)  a - — 2k(|~z°-)  + image  term.  (4-5) 

o t G q K 

V2<f  at  a later  time  may  also  be  predicted.  With  the  new  V2$  distri- 
bution we  can  find  a new  value  of  3/3t(V24>),  and  so  advance  for  <}>, 
the  solution  in  time.  This  "finite  difference  method"  was  done  by 
Illingworth  (1972). 

Estimates  for  time  constant  of  recoveries  are  thus  20.4  sec  at  18km, 
8.4  sec  at  24  km  and  4.6  sec  at  35  km  for  a discharge  at  6 km  regardless 
of  the  amount  of  the  discharge.  This  explains  the  more  rapid  recoveries 
at  higher  altitudes  observed  above  clouds. 

Recent  observations  of  balloon  flights  by  Winn  and  Byerly  (1975) 
also  indicated  faster  recoveries  of  electric  fields  at  ground  than  at 


k 
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balloon  altitude.  Furthermore,  they  described  the  pattern  of  field 
change  as  a nearly  linear  function  of  time.  This  is  not  a direct  impli- 
cation of  (4-5).  However,  there  are  also  other  processes  affecting  the 
field  recoveries  in  addition  to  the  regeneration  of  the  charge  destroyed. 
Examples  are  the  rearrangement  of  space  charge  in  the  conductivity 
gradient  between  the  cloud  and  the  air,  the  effects  of  local  point 
discharges  and  reversal  distance  complications  for  cloud  discharges. 

These  may  all  change  the  predicted  shape  of  the  recovery  curve. 

Furthermore,  we  observed  a definite  correlation  between  the  amount 
of  discharge  and  recovery  time  constant  (Figure  4-14) . The  time  con- 
stant for  field  recovery  is  proportional  to  the  amount  of  discharge. 

The  simple  linear  correlation  can  be  written  as 

AE  = 485.4  + 25.7  x t (4-6) 

max 

where  AE  is  in  mV/m  and  T in  seconds.  The  correlation  coefficient 
max 

r (a  measure  that  shows  how  well  the  linear  curve  actually  fits  the  data; 
r = 1 means  perfect  correlation  and  r = 0 is  totally  uncorrelated) 
is 

m v 2 
ax 

r = r 

v i 
ay 

where  m is  the  slope,  v 2 is  the  variance  in  the  x values,  and 

ax 

v 2 is  the  variance  in  the  y values.  In  our  case 
ay 

25.7  x v 
ar 

VaAE 
= 0.83. 
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The  lack  of  knowledge  of  the  condition  inside  a storm  prevents  us  from 
pursuing  explanation  for  this  observation  further. 

A better  and  more  thorough  system  of  solving  the  problem  is  defi- 
nitely necessary.  Detail  observation  of  lightning  field  recoveries  at 
different  altitudes  will  also  help  in  solving  the  mystery  of  charge 
mechanism  during  the  active  period  of  a thunderstorm. 


ft— 


Figure  4-1 

Map  of  Canada  giving  locations  of  the  balloon 
launch  sites.  (Moser  and  Manka,  1971) 
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Figure  4-3a 


A portion  of  raw  data  indicating  the  presence  of  thunderstorm. 


Figure  4- 3b 

Another  portion  of  raw  data  indicating  the  presence  of 
lightning  activity 
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VERTICAL  ELECTRIC  FIELD  EXPECTED 
AT  BALLOON  ALTITUDE 


Figure 


EXPECTED  ELECTRIC  FIELD  IN  THE  HORIZONTAL  PLANE 
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Figure  4-8 

Vector  electric  fields  during  the  thunderstorm  from 
0500  to  0630  UT  observed  by  the  balloon  sensors  at 
90,000  ft. 
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Figure  4-9 

Observed  vertical  electric  field  as  compared  to  different 
sources: 

(1)  field  due  to  thundercells  (upward), 

(2)  shielding  field  (downward),  and 

(3)  normal  background  field  (downward). 

The  addition  of  the  three  sources  shows  agreement  with 


the  observed  field. 


Figure  4-11 

(a)  Electric  field  components  vs.  time  for 
a return  stroke. 


(b)  Total  electric  field  intensity  vs.  tim 
at  several  distances  due  to  the  stroke 
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Figure  4-13 


Changes  in  the  vertical  fields  indicating  the  presence 

of  lightning  flashes.  A E is  the  maximum  change 

max 

in  field.  T is  the  relaxation  time  for  field  recovery. 


linear  corelation  between  E and 

max 
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CHAPTER  5 
DISCUSSION 

I.  An  Outline  of  Electrical  Structures  In  the  Atmosphere,  Ionosphere, 
and  Magnetosphere 

In  the  previous  chapters,  we  have  studied  the  property  of  air  con- 
ductivity and  discussed  how  the  atmospheric  electric  fields  are  mapped  at 
different  altitudes.  Comparisons  between  the  model  and  balloon  data  con- 
firmed the  relative  position  between  the  balloon  and  the  thunderstorm. 

The  charge  concentration  inside  different  cells  was  thus  calculated. 

Furthermore,  the  mapping  technique  was  also  extended  into  the  ionosphere. 

In  this  chapter,  we  will  review  the  large  scale  coupling  problem  in  both 
theoretical  and  experimental  aspects  and  discuss  the  role  of  mapping  tech- 
nique in  solving  the  problem. 

Since  Benjamin  Franklin  discovered  electricity  in  lightning,  the 
study  of  atmospheric  electricity  has  been  concerned  mainly  with  the  fol- 
lowing problems: 

(1)  the  nature  and  origin  of  the  atmospheric  conductivity, 

(2)  the  origin  and  maintenance  of  earth's  negative  charge,  and 

(3)  the  source  and  development  of  thunderstorm  charges. 

Meanwhile,  the  discovery  of  the  regular  diurnal  variations  in  the  geo- 
magnetic field  has  made  the  subject  of  space  electricity  also  the  con- 
cern of  ionospheric  and  magnetospheric  researchers.  It  was  the 
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investigations  of  whistler  phenomena  that  led  to  the  hypothesis  that  there 
exists  a coupling  relationship  between  the  atmosphere,  the  ionosphere  and 
the  magnetosphere.  General  electrodynamic  mechanism  for  each  region  is 
illustrated  for  an  idealized  height  distribution  in  Figure  5-1  (Obayashi 
and  Maeda,  1965). 

In  the  atmosphere,  there  exists  a downward  electric  field  during 
fair  weather  conditions  since  the  earth  is  negatively  charged.  The  elec- 
tric field  is  maintained  mostly  by  convection  currents  provided  by  source 
currents  from  storm  weather.  The  normal  field  intensity  is  about  130  V/m 
near  ground  but  decreases  with  increasing  altitude  due  to  the  increase  of 


1 
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conductivity  at  higher  altitudes.  The  total  potential  difference  between 
the  earth  and  the  ionosphere  is  in  the  order  of  4 x 105  volts.  The  cur- 
rent-electric field  relationship  in  this  region  can  be  written  as 


(5-1) 


At  ionospheric  heights,  the  electric  field  is  generated  by  con- 
vective motions.  The  sources  at  different  latitude  regions  are: 

(1)  At  lower  and  middle  latitude  regions,  the  source  is  the 
interaction  of  tidal  winds  with  geomagnetic  field  and 
charged  particles  in  the  upper  atmosphere,  the  so-called 
"dynamo  theory"  (Chapman  and  Bartels,  1940). 

(2)  At  higher  latitude  regions  (aurora  and  polar  regions),  the 
plasma  flow  and  electric  field  pattern  is  shown  in  Figure 
5-2.  The  two  possible  sources  for  the  antisunward  plasma 
flow  over  the  polar  cap  and  the  sunward  return  flow  in  the 
aurora  zone  are:  (a)  viscous  interaction  (Axford  and  Hines, 
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1961),  and  (b)  magnetic  merging  (Dungey,  1961). 

The  average  magnitude  of  the  field  varies  between  0 and  20  mV/m, 
depending  on  the  conductivity  in  the  region,  the  latitude  location  and 
local  time.  Large  diurnal  changes  of  amplitude  and  direction  are  the 


general  rule  in  the  ionosphere.  The  electric  field  which  drives  the 
current  consists  of  two  distinct  parts:  one  due  to  the  convection  motion, 


producing  an  induction  field  E^  = v x B,  and  the  other  is  the  electro- 
static field  Eg  resulting  from  the  accumulation  of  polarization  charges. 
The  expression  can  be  written  as 


j • (a)  = E + v x B 


Es  = j • (g>  - V X B 


(5-2) 


In  the  magnetosphere  the  electric  fields  are  generated  by  corota- 
tional  motions  within  the  plasmasphere  and  convective  motions  in  the 
plasmasheet  and  high  latitude  lobes.  The  electrostatic  field  is  such 
that  it  compensates  the  Lorentz  force  in  order  to  preserve  the  convective 
motion.  The  representative  magnitude  is  in  the  order  of  1 mV/m.  The 
corresponding  expression  in  this  region  is 


j ’ (o)"1  = 0 


E = - v x B 
s 


(5-3) 
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II . Techniques  of  Electric  Field  Mapping 

(i)  Mapping  between  the  Magnetosphere  and  Ionosphere 

Wolf  (1975)  has  given  a thorough  review  on  the  coupling  problem  of 
the  ionosphere  and  magnetosphere.  Magnetic  field  lines  in  the  higher 
ionosphere  and  magnetosphere  can  be  considered  to  have  high  enough  con- 
ductivity, at  most  times  and  places,  such  that  a reasonable  zero  order 

picture  of  circulation  in  the  ionosphere  and  magnetosphere  can  be  obtained 

•>  ->• 

with  the  assumption  that  E • B = 0.  In  the  region  of  dipolar  magnetic 
field  lines,  Mozer  (1970)  developed  the  relationships  of  electric  fields 
between  the  two  regions  as 


gionosphere 

east 

^.equat.  plane 
east 


(5-4) 


E 


ionosphere 

radial 

equat.  plane 
radial 


(5-5) 


where  the  L-value  of  a particular  point  is  the  distance  from  the  center 
of  the  earth  at  which  the  field  line  through  that  point  crosses  the 
equatorial  plane.  The  unit  of  L is  in  earth  radius. 

Mapping  of  fields  in  the  outer  magnetosphere  is  a more  complicated 
problem  and  will  not  be  discussed  here  (for  more  references,  see  the 
review  paper  by  Wolf,  1975). 

Components  of  electric  field  parallel  to  magnetic  field  lines  are 


the  least  understood  even  though  their  existence  has  been  suggested  by 
data  from  rocket-borne  double  probe  experiments  (Mozer  and  Fahleson,  1970; 
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Kelley  et  al.,  1971).  More  experiments  must  be  done  in  order  to  produce 
information  on  the  average  global  distribution  of  total  field  aligned  potentia 
drops  and  thus  to  determine  whether  these  fields  play  an  important  role 
in  large  scale  ionosphere-magnetosphere  dynamics. 

(ii)  Mapping  between  the  Ionosphere  and  Atmosphere 

Unlike  the  two-way  relationship  of  electric  fields  between  the 
ionosphere  and  the  magnetosphere  as  expressed  in  (5-4)  and  (5-5) , the 
field  mapping  from  ionosphere  down  to  the  atmosphere  is  quite  different 
from  that  of  the  atmosphere  into  the  ionosphere.  The  main  reasons  are: 

(1)  In  the  atmosphere,  the  weather-induced  fields  are  due  to 

charge  sources.  In  the  ionosphere,  electric  fields  are 

*> 

mostly  induced  from  v x B convections.  Because  of  the 
high  conductivity,  regions  of  net  charges  cannot  exist 
in  the  ionosphere. 

(2)  In  the  atmosphere,  the  conductivity  of  air  is  isotropic 
and  can  be  expressed  as  a scalar.  In  the  ionosphere, 

the  conductivity  becomes  anisotropic  and  the  presentation 
is  in  tensor  form. 

Recent  work  on  downward  field  mapping  from  the  ionosphere  into 
the  atmosphere  include  Mozer  (1971)  ano  Volland  (1972).  Assuming  scalar 
conductivity  up  to  100  km  where  ionospheric  fields  originate,  both  of 
them  showed  that  these  fields  can  map  into  the  lower  atmosphere  (^30  km) 
without  much  attenuation.  Using  numerical  techniques,  Park  (to  be  pub- 
lished) gave  a more  thorough  treatment  taking  into  consideration  of  the 
tensor  presentation,  possible  horizontal  gradient  and  other  irregularities 
of  ionospheric  conductivity.  Thus  more  generalized  results  are  obtained. 
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Upward  mapping  of  atmospheric  electric  fields  into  the  ionosphere 
was  reviewed  in  Chapter  2.  In  this  thesis  we  developed  an  analytic  solu- 
tion for  mapping  of  thunderstorm  electric  fields  to  different  altitudes, 
including  balloon  altitudes  and  ionosphere.  The  results  at  balloon  alti- 
tude agree  with  measured  values.  Near  the  boundary  between  the  atmosphere 
and  the  ionosphere  ( ^ 70  km),  the  thunderstorm  induced  fields  are  in  the 
same  order  of  magnitude  as  ionospheric  fields  — several  to  tens  of  milli- 
volts per  meter.  In  the  ionosphere,  the  fields  attenuate  much  faster  so 
that  thunderstorm  induced  fields  will  fall  below  the  instrument  detection 
range  at  ^ 100  km. 

It  seems  paradoxical  that  ionospheric  fields  can  map  down  to  the 
atmosphere  but  the  atmospheric  fields  dissipate  rather  quickly  in  the 
ionosphere.  The  explanation  for  this  is  simply  that  ionospheric  conduc- 
tivity is  much  greater  than  atmospheric  conductivity;  the  electric  field 
is  mapped  more  efficiently  along  a path  of  decreasing  conductivity  (k < 0) 
than  vice  versa. 

Ill . Experimental  Techniques 

The  advantage  of  vector  electric  field  measurements  above  thunder- 
storms is  they  avoid  all  ambiguities  near  ground;  for  example,  corona 
discharges  near  trees  and  other  ground  conductors  and  the  anomalies  of 
conductivity  between  thunderclouds  and  ground.  As  a result,  measurements 
above  thunderstorms  can  give  us  first-hand  information  on  the  time  history 
and  electrical  structure  of  the  storm  and  thus  help  in  monitoring  the 


electrification  processes  inside  thunderclouds. 
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Another  advantage  of  field  measurements  at  altitudes  between  30 
and  40  km  is  that  electric  fields  of  both  atmospheric  and  ionospheric 
origins  can  be  observed  and  studied.  They  can  then  provide  valuable 
evidence  in  the  ionosphere-atmosphere  coupling  problem. 

Popular  vehicles  for  electric  field  measurements  at  higher  alti- 
tudes include  (1)  airplanes,  (2)  barium  cloud  releases,  (3)  rockets  and 
(4)  balloons.  Airplanes  provide  the  best  control  of  detecting  path. 

Their  major  problem  lies  in  the  determination  of  the  enhancement  factors 
since  the  payload  is  usually  mounted  at  the  tip  of  the  plane  where  field 
distortion  is  the  worst.  Barium  cloud  releases  are  easy  to  perform,  yet 
its  application  is  mainly  on  the  v x B,  convection,  electric  fields. 

Storm  induced  fields,  since  they  do  not  always  cross  the  magnetic  field 
lines,  cannot  be  accurately  determined.  Rockets  are  capable  of  accurate 
measurements.  However,  their  flight  durations  are  obviously  not  suffi- 
cient for  studying  the  time  variation  of  the  fields.  Balloons,  on  the 
other  hand,  even  though  it  is  not  possible  to  control  their  exact  path, 
compensate  for  most  of  the  disadvantages  of  the  other  vehicles  and  thus 
are  good  payload  carriers.  Of  course,  in  order  to  get  a complete  picture 
of  the  electrical  structure  of  upper  atmosphere,  simultaneous  measurements 
from  satellite  and  ground  bases  are  also  necessary  in  addition  to  balloon 
measurements. 

IV . Future  Prospects 

In  short,  we  have  shown  the  mapping  technique  from  the  atmosphere 
into  the  ionosphere.  Our  conclusion  is  that  thunderstorm  induced  elec- 
tric fields,  after  mapping  into  the  ionosphere,  are  much  less  than  field 
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produced  by  ionospheric  or  raagnetospheric  processes.  Nonetheless,  such 
fields  might  be  responsible  for  formation  of  localized  electron  density 
enhancements,  the  so-called  "magnetospheric  ducts"  (Mozer,  1971). 

Another  aspect  requiring  research  is  the  VLF  range  of  the  lightning 
induced  fields;  even  though  electrostatic  fields  dissipate  fast  in  the 
ionosphere,  the  radiation  field  is  still  a possible  excitation  source  of 
geomagnetic  pulsations  in  the  magnetosphere  (suggested  by  Dejnakarintra 
and  Park,  1974). 

Furthermore,  the  subject  of  solar-terrestrial  effects  on  the 
weather  has  drawn  more  and  more  interest  and  attention  these  days  (Mark- 
son,  1971,  1975;  Roberts  and  Olson,  1973;  Silverman  and  Korff,  1975). 
Mapping  techniques  between  magnetosphere,  ionosphere  and  atmosphere  will 
certainly  be  a first  necessary  step  in  setting  up  a theoretical  model. 

It  is  certainly  encouraging  that  the  science  of  atmospheric  elec- 
tricity is  now  experiencing  more  interactions  with  other  disciplines  of 
space  studies.  It  will  work  to  the  benefit  of  all  of  science  that  these 
disciplines  are  now  cooperating  after  so  many  years  of  isolated  research. 
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Mathametical  Appendix 


An  Alternate  Method  of  Solving 
the  Potential  Problem  in  the  Atmosphere  with  a Point  Source 
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The  Laplace  equation  to  be  solved  is 


V2  + 2k  -^-1 4>  = 0 


with  boundary  condition  of  charge  at  z^. 

In  the  cylindrical  coordinates: 


Assuming  azimuthal 
gives 


1 _3_ 
r 3r 


- If] 


symmetry,  separation  by  means  of  the  product  function 


i 

f 

We  obtain 

i 


i 


<J>(r,z)  = R(r)Z(z) . 


d2R  1 dR 
dr2  rdr+ 


i2R  = 0 


i 


d2Z 

dz2 


+ 


2k  - m2Z  = 0 
dz 


(A-l) 

(A-2) 


where  m is  a separation  parameter.  Equation  (A-l)  is  known  as  Bessel’s 
equation. 

The  solution  of  (A-l)  is  of  the  form 


R(r)  = C,J0(mr)  + C2Y0(mr) 

where  J0(  > and  Y0(  ) are  Bessel  functions  of  first  and  second  kind 
with  zero  order.  Boundary  condition  here  requires  that  as  r -*•  0,  R(r) 
remains  finite 


C - 0 


such  that 


R(r)  « Cj  J0  (n>r) 

The  solution  of  (A-2)  is  of  the  form 


~ f U ~ 

Z(z)  = C3exp  |/m2+k2  - k z-z0 


C4expj^  |-/mz-Hnz  - kj  [z"z0] 


(A- 3 


Boundary  conditions  in  this  case  require  that 

as  z -*■  Z(z)  0 such  that  C3  = 0 for  z > zQ 
and  as  z -*•  -•»,  Z(z)  ■*  0 such  that  C4  = 0 for  z < zQ 


Z(z) 


C4exp 

/m2+tn2  - kj 

z-z 

0 

] 

C3exp 

/m2+kz  - k 

_ 

z"2o 

J 

] 

for 


for 


z > zQ 


z < z, 
0 


(A-4 


(A-5 


where  z > zQ  and  z < z0  refer  to  above  and  below  the  charge  source 
respectively.  (We  are  interested  in  the  z > zQ  case.) 

For  z > z0,  we  can  write 

<A-6] 


y A J (mr)  exp  - /mz+k2  - k 
m 0 


z-z„ 


and  for  z < z^, 


= I AmJfl(mr)  expp  [ /m2+k2  - k z-zQj  ~]  . (A-7} 
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In  order  to  find  Affl,  we  apply  boundary  conditions  at 
of 


Z = Z„  by  means 


fOO  rao 

5(r)  <5(z-z.)  2tvr 

J0  J0 


dr  dz  * 1 


or 


fu 

Jn 


r ) 2trr  dr  = 1 

as  shown  in  Phillips  and  Panof sky  (1962) . 

The  flux  condition  at  the  plane  z = z is 


(A-8) 


-Jz=z„  o J 


p dV 


which  implies  that 


2 I &&  AmJfl(mr)  = 
m=o 


ra 


(A-9) 


Multiplying  both  sides  of  (A-9)  by  rJJm’r)  and  integrating  over  r <0~) : 

f 


by  using  the  relationship  (A-8). 
Since  Jo(0)  - 1,  and 


l X Vkx)  Jp(k'x)  dx  = i 6(k'-  k) 


or 


f x (J  <*«)) 

Jq  P 


2 a 1 

dx  * r . 


1 
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we  obtain 


or 


2/m2+kJ  • A • — 


m m 2tre  ( 


m “ 4Tre0  * /m2+k2 


Now  for  z > z. 


. . , -k(z~20)  r 

4>(r,z)  = -a- 


4ire, 


fCO 

m j , x -(z-z0)/m7HhkT 
J <mr)  e 


dm 


By  using  the  relationship  in  Gradshteyn  and  Ryzhik  (1955,  6.616) 

| J0(r/m'  - 1 ) exp(-zm')  dm' 


exp  f-/r2-*-zz  ) 
/r2+z2 


we  get 


, , . -k/r2+(z-z  ) 2 

*fr  -■>  - 9 .-k(z-zo)  e 0 _ 

*(r’z)  47ie0  e /r»  + (z  - z0Y2  * 

Introducing  an  image  charge  of  +q  at  -z„  to  match  the  boundary 
condition  on  the  ground  such  that  <J>  = 0 at  z = 0, 


-k/r4+(z-zn ) 2 -k^r2+(z+zn) 2 


/r*  + (z  - z0)2 


/r2  + (z  + zo ) 2 
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